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ABSTRACT 


CEINCI-LAB is a computer tool to perform seismic and static calculation of 
structures, it is oriented to academia and research because it allows the user to 
understand the data entry, the calculation sequence and the results of the analysis, 
also It can be applied in consultancies freely. 


This article presents new functions of CEINCI-LAB oriented to the analysis 
and design of steel frames; in the seismic analysis, displacements, drifts and shear 
forces are calculated on each floor. In the static analysis, the gravitational and 
seismic loads are applied to obtain the forces in the elements, the load 
combinations of the Ecuadorian Construction Standard (NEC-15) are applied and 
the demand envelopes are obtained to compare them with the capacities of the 
elements. It is also determined that the beam / column ratio is less than unity to 
ensure adequate seismic behavior with the dissipation of energy in the beams. 


These functions allow the user to easily view the results in graphs with 
colors depending on the demand versus capacity relationship of the elements, in a 
similar way to that presented in commercial programs such as ETABS or 
SAP2000, in addition the user can adapt the functions to particular needs thanks to 
the goodness of free code. 


keywords: CEINCI-LAB computer system. Seismic Analysis. Static Analysis. Steel 
structures. NEC - 15. 


2 Brian Cagua, Roberto Aguiar y Julia Pilatasig 


NUEVAS FUNCIONES DE CEINCI-LAB PARA EL 
ANALISIS Y DISENO DE PORTICOS DE ACERO 
ACORDE A LA NEC-15 


RESUMEN 


El Sistema de Computación CEINCI-LAB es una herramienta informática 
para realizar cálculo sísmico y estático de estructuras, se orienta a la academia e 
investigación debido a que permite al usuario entender el ingreso de datos, la 
secuencia de cálculo y los resultados del análisis, también se puede aplicar en 
consultorías de manera libre. 


En este artículo se presenta nuevas funciones de CEINCI-LAB orientadas 
al análisis y diseño de pórticos de acero; en el análisis sísmico se calcula 
desplazamientos, derivas y fuerzas de corte en cada piso. En el análisis estático 
se aplican las cargas gravitacionales y sísmicas para obtener las fuerzas en los 
elementos, se aplica las combinaciones de cargas de la Norma Ecuatoriana de la 
Construcción (NEC-15) y se obtiene las envolventes de la demanda para 
compararlas con las capacidades de los elementos. También se verifica que la 
relación viga/columna sea menor a la unidad para que se asegure un 
comportamiento sismorresistente adecuado con la disipación de energía en las 
vigas. 


Estas funciones permiten que el usuario pueda visualizar los resultados de 
manera sencilla en gráficos con colores en función de la relación demanda versus 
capacidad de los elementos, de manera similar a la que se presenta en programas 
comerciales como ETABS o SAP2000, además el usuario puede adaptar las 
funciones a necesidades particulares gracias a la bondad que presenta el código 
libre. 


Palabras Claves: Sistema de computación CEINCI-LAB. Análisis Sísmico. 
Análisis Estático. Estructuras de Acero. NEC — 15. 


1. INTRODUCCIÓN 


Importantes contribuciones han sido desarrolladas por los investigadores 
del blog de CEINCI-LAB para facilitar el análisis sísmico de estructuras en tres 
dimensiones (Aguiar et al. 2019); para la generación automática y análisis de 
estructuras de hormigón armado orientado al cálculo de la vulnerabilidad sísmica 
de una gran cantidad de estructuras (Pilatasig et al. 2020). 


También se ha incursionado y se continúa trabajando en el acoplamiento 
entre los sistemas de computación CEINCI-LAB y OpenSees. Este acoplamiento 
permite analizar estructuras en forma amigable y sencilla. Una de las fortalezas de 
CEINCI-LAB es la generación de estructuras y la presentación gráfica de 
resultados de una manera muy fácil y la fortaleza de OpenSees son las funciones 
con las cuales se realiza análisis no lineal, fundamentalmente. Cagua et al. 2020; 
Aguiar et al. 2020. 
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Ahora bien, en el Ecuador y en algunos países que tienen una alta 
peligrosidad sísmica se están construyendo una gran cantidad de edificios en 
acero, compuestos por pórticos sin arriostramiento, como los indicados en las 
fotografías de la figura 1; la unión de las columnas y vigas de las estructuras 
presentadas son soldadas. Todavía no está generalizado el uso de conexiones 
precalificadas, con las que se tiene conexión de momento. 
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Figura 1 Estructuras de acero en el Valle de los Chillos; a la izquierda un centro 
comercial y a la derecha una edificación de apartamentos. 


Este artículo está orientado al análisis y diseño de estructuras de acero 
conformadas por vigas y columnas con conexiones que tienen capacidad de 
momento, para el efecto se han implementado nuevas funciones al sistema de 
computación CEINCI-LAB. 
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Figura 2 Filosofia de diseno en estructuras de acero, sin sistemas de control 
activo como aisladores sismicos. Adaptado de: Guia de diseno 3 NEC. 





Vale la pena recordar la filosofia de diseno de estructuras de acero sin 
sistemas de control (aisladores sismicos o disipadores de energia). Para el efecto, 
en la grafica de la derecha de la figura 2, se muestra la capacidad sismica 
resistente, que relaciona el cortante basal con el desplazamiento lateral en el tope 
del edificio, de dos estructuras, identificadas como 1 y 2. 
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La Estructura 1 de la figura 2, ha sido diseñada utilizando un espectro 
elástico, de tal manera que ante un sismo determinado la demanda sísmica, 
desarrolla un cortante basal V, asociado a un desplazamiento lateral máximo Ag. 
La rigidez de la estructura es la pendiente de la recta y se ve que esta es 
constante, de tal forma que la estructura no sufre ningún daño, pero al ser 
diseñada para un valor V, muy alto las dimensiones de sus elementos son grandes 
y por ende se trata de una edificación costosa y a lo mejor no es tan funcional. 


Por otro lado, la estructura 2 (figura 2), fue diseñada para un espectro que 
se obtiene dividiendo el espectro elástico para el factor de reducción R con el cual 
se tiene un espectro inelástico con ordenadas espectrales bajas que conduce a 
tener elementos estructurales de menores dimensiones. Esto implica que durante 
el mismo sismo la estructura va a disipar energía por agotamiento de elementos a 
flexión debido a las conexiones que se habló anteriormente. 


Se realiza el análisis sísmico utilizando el espectro inelástico, con lo que 
se hallan las fuerzas horizontales F; que se muestran en el Pórtico de la figura 2, la 
suma de estas fuerzas laterales, dan origen al cortante basal V,. En la figura de la 
derecha se observa el cortante basal V, asociado a un desplazamiento lateral A,. 
Este punto es el límite del rango elástico de la Estructura 2. Para encontrar el 
desplazamiento inelástico A, que realmente va a tener la estructura durante el 
sismo severo se multiplica el desplazamiento elástico A, por el factor C4 que 
recomienda el ASCE para diferentes tipologías estructurales. 


Con esta introducción se debe indicar que la filosofía de diseño en acero 
comprende los denominados Sistemas Resistentes a cargas Sismicas (SRCS), 
que incluyen a todos los elementos de acero estructural y conexiones de un 
edificio que han sido diseñadas para resistir cargas laterales y verticales, y los 
elementos que poseen la característica de resistir solo cargas verticales. (Esta 
forma de diseño también se aplica a estructuras de hormigón armado y se las ve 
un poco diferentes, un ejemplo se tiene en las estructuras de la Universidad de 
Fuerzas Armadas ESPE, en Santo Domingo). 


El procedimiento específico de diseño de estructuras de acero depende de 
las particularidades de cada proyecto, pero todas las estructuras deberán ser 
diseñadas bajo el mismo procedimiento general (SGR, MIDUVI, 2016), como se 
muestra a continuación en la figura 3. 


Aparentemente es bastante sencillo determinar las cargas gravitacionales 
y sísmicas, pero en realidad no lo es, por eso se recomienda leer las 
especificaciones que se indican en la parte superior de la figura 3. 


En este artículo se presenta el análisis sísmico en dos dimensiones, para 
poder ilustrar el análisis y diseño sísmico de una estructura de acero; con pocos 
cambios se puede adaptar al análisis en tres dimensiones. Se recomienda para 
ello en una primera etapa hacer el cálculo espacial de estructuras regulares con un 
modelo de un grado de libertad por piso en el centro de masas e incluir en el 
estudio la torsión accidental y después pasar al análisis sísmico en tres 
dimensiones. 
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Figura 3 Procedimiento general para el diseño de estructuras de acero, según 
NEC 2015 


Existen varios sistemas resistentes de cargas laterales utilizados en la 
construcción con acero estructural, a continuación, se describen los tres sistemas 
más comunes. 


Pórticos Resistentes a Momento 


La característica principal de este tipo de pórticos es la formación de 
rótulas plásticas de manera secuencial en las vigas, como último mecanismo de 
disipación de energía se permite la formación de rotulas plásticas en la base de las 
columnas de planta baja como se muestra en la figura 4. Además, en estos 
pórticos, las conexiones entre vigas y columnas son rígidas. 


NUDO PLASTICO e NUDO PLASTICO 





NUDO PLÁSTICO 


NUDO PLÁSTICO R’ E ra NUDO PLÁSTICO 


(a) (b) 


Figura 4 a) Estructura sin deformar con conexión rígida; b) Estructura deformada. 
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A continuación, en la figura 5, se describe algunos de los parámetros que 
se deben tomar en cuenta para asegurar el comportamiento sismorresistente 
adecuado de los pórticos resistentes a momento. 






Escoger los lugares en donde se formarán las rótulas plásticas. 


Diseñar y detallar los lugares escogidos para que puedan soportar 
deformaciones considerables y no ocurra un fallo por fractura del 
componente o por inestabilidad. 






Diseñar el resto de los elementos del pórtico para que sean más resistentes 


que las regiones en donde se formarán las rotulas 





Figura 5 Parámetros de un diseño sismorresistente. 


Se puede asegurar la formación de la rótula plástica en el lugar “escogido”, 
mediante un diseño de la conexión denominado RBS, esto implica secciones de 
patín reducido en las vigas, es decir, que se puede acortar el patín con un radio de 
dimensiones controladas y de esta manera reducir la capacidad a flexión del 
elemento en esa zona. 


Una de las desventajas que se puede mencionar, es que en este tipo de 
sistemas la rigidez lateral es relativamente baja, lo que se traduce como derivas 
mayores a las de otros sistemas. (SGR, MIDUVI, 2016). 


En este artículo se realiza el análisis y diseño estructural de un pórtico de 
acero empleando el Sistema de Computación CEINCI-LAB, se detallan las nuevas 
funciones desarrolladas para el análisis sísmico, específicamente para calcular 
desplazamientos, derivas y fuerzas de corte en cada piso mediante el método 
Estático Equivalente y Modal Espectral. Además, para el análisis estático se 
aplican las cargas gravitacionales y sísmicas obteniendo las fuerzas internas en 
los elementos, se consideran las combinaciones de cargas de la Norma 
Ecuatoriana de la Construcción (NEC-15) y se obtiene las envolventes de la 
demanda. La demanda se compara con la capacidad (a carga axial, corte y 
momento) de los elementos y se calcula la relación viga/columna, se debe cumplir 
con un valor menor a 1 para que se asegure un comportamiento sismorresistente 
adecuado con la disipación de energía en las vigas (columna fuerte/viga débil) 
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2. GEOMETRÍA Y CONSIDERACIONES DE CARGAS 


De manera acádemica se realiza el análisis de un pórtico plano, esto 
debido a que es más sencillo detallar los calculos y también permite que el lector 
asimile la información de mejor forma. La geometria del pórtico considerando 3 
vanos y 4 pisos se muestra en la figura 6. Se considera que las cargas 
gravitacionales se aplican a las vigas del pórtico mediante una distribución 
uniforme por razones de simplificación en el cálculo. 
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Figura 6 Geometria del portico analizado 


Para el diseno se considera las cargas gravitacionales sobre el sistema de 
piso; la carga muerta se compone del peso propio de la losa, vigas y columnas, 
ademas de la carga permanente, para los entrepisos se considera de 500 kg/m? y 
para cubierta de 350 kg/m. La carga viva depende de la ocupación, en este caso 
se asigna 250 kg/m” a los entrepiso y 100 kg/m? a cubiertas (tabla 1). 


Tabla 1 Cargas gravitacionales consideradas 


Parámetro | Valor | Unidad | Observación 
Peso de placa = | 637 | wm | O O 
[Volumen hormigón == [007854 | mm [SSCS 
[Peso de hormigón == [19089 | m [SSCS 
[Otros (Instalaciones) = | 6 | ome | SSCS 
Enlucido masilado= | a | kom | 2omdeenuido 


Asumido (Depende de proyecto 
arquitectónico) 


Peso mampostería = kg/m? 


P i | 
eso propio de elementos kg/m? Asumido 
estructurales = 


Carga permanente cubieria= | 150 | kom [SSCS 
[Carga viva pisos | 250 | kgm | —NEC_SECG2015, pig. 25 | 
[Carga viva cubieraz | 100 | kgm? | NEC SE 0G 2015 pág. 25 | 
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3. COMPORTAMIENTO DE LOS ELEMENTOS DE ACERO 
SEGUN SU ESBELTEZ 


Las vigas y columnas de un pórtico de acero se encuentran sometidas a 
efectos de flexión, la capacidad de estos elementos de soportar estas 
solicitaciones es en función de la compacidad de las secciones, como se indica en 
la figura 7, es decir de la clasificación de las secciones en función de sus 
relaciones de ancho versus espesor. 
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Figura 7 Capacidad a flexión de una sección en función de su compacidad 


= Secciones con alta ductilidad (1). Se puede alcanzar un momento 
resistente igual al momento plástico, además podría desarrollar grandes 
deformaciones unitarias sin falla prematura por inestabilidad local 
(pandeo). La relación de anchura-espesor À < Apa- 

= Sección moderadamente dúctil (2). Sección que tiene poca probabilidad 
de presentar pandeo local, no obstante, puede mantener deformaciones 
unitarias tan altas como en el caso anterior. En este caso: Apa <A < Àp. 

= Sección no compacta (3). Pasado del límite de fluencia se pierde su 
capacidad de carga a causa del pandeo local (falla inelástica). A > 2, 

= Sección esbelta (4). La sección no desarrolla su máxima capacidad de 
carga debido a una falla por pandeo, esto limita también su capacidad de 
deformación. En este caso 1 > 1, 


El pandeo es un fenómeno conocido como la inestabilidad elástica que 
puede darse en elementos comprimidos, en especial cuando sus elementos son 
esbeltos, y que se manifiesta por la aparición de desplazamientos importantes 
transversales a la dirección principal de compresión (figura 8). 


INICIO DEL VOLADO 


FIN DEL VOLADO 





Figura 8 Pandeo lateral torsional en un volado. 
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Los efectos de pandeo o inestabilidad podrían presentarse de manera 
local en una sección, como se ilustra en la figura 9. 


FALLA POR PANDEO LOCAL 





TENSIÓN ! 


a) b) 
Figura 9 a) Pandeo local del patin b) Pandeo local del alma. 


La función Compacidad_Elementos determina la relación de aspecto 
b/t; (ala) y h/t,, (alma); los valores límites que se programaron son los del AISC 
341-16. El programa clasifica a los elementos en: HD altamente dúctil; MD 
moderadamente dúctil y ND poco dúctil. Más adelante se habla sobre esta función. 


3.1 CAPACIDAD A FLEXIÓN DE ELEMENTOS DE ACERO 


Al aplicar un momento flector a una viga, la sección transversal de la 
misma se deforma y debe resistir esfuerzos internos debido a esta solicitación. En 
primera instancia se supone que el material (acero) se comporta en su rango 
elástico y las deformaciones que se observan en la sección son elásticas, es decir 
que al retirar la carga, la sección regresa a su condición inicial sin sufrir daños; sin 
embargo, si el momento aplicado (solicitación) supera el esfuerzo de fluencia, es 
decir si la sección supera el rango elástico, la viga contempla deformaciones 
plásticas que al retirar la carga ya no se recuperan, en consecuencia la sección 
queda “deformada”. Una secuencia de los esfuerzos en función del momento 
aplicado se tiene en la figura 10. 
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M Elastic M 
> 100 ) M 
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6, x E, g ~ F G, 7 F. 
M<M, M,< MSM, M= M, 


Figura 10 Esfuerzos internos en una sección en función del momento aplicado. 


Las vigas de acero pueden llegar a tener el comportamiento mostrado en 
la figura 11; es importante evitar la inestabilidad global y local de las vigas para 
poder desarrollar su capacidad de plastificación (generalmente se recomienda que 
las secciones sean compactas y que los elementos tengan soporte lateral para 
tener este comportamiento, pero ante un terremoto es preferible que trabajen en el 
rango elástico). 
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Figura 11 Diagrama momento versus curvatura en una sección l. 


= Diseño elástico: 


El momento elástico (Me) es el momento para el cual la fibra más alejada 
del eje neutro alcanza el esfuerzo de fluencia (Fy), este valor se encuentra 
relacionado con el módulo de la sección elástico, que generalmente está detallado 
en catálogos comerciales de perfiles y presenta la nomenclatura “W” o “S”. El 
módulo de sección es la relación de la inercia (l) para la distancia a la fibra más 
externa del perfil desde al eje neutro, en perfiles simétricos respecto al eje “x” esta 
distancia es igual a la del centro de gravedad (ycg) como se muestra en la figura 
12. 


AY 





máx Fy 
I 
xx 
d úx Wy —x ~~ Saa = 
Ycg 
eje neutro — 
Po > Mey-x = Fy * Sy_x 


Figura 12 Esfuerzos internos en una sección sometida al momento elástico. 


= Diseño plástico: 


Una sección para poder desarrollar la plastificación debe cumplir con 
ciertos criterios para garantizar su respuesta inelastica, estos criterios son: 


= Calidad del acero utilizado para poder alcanzar ductilidad. 

= Relaciones geométricas de su sección para evitar el pandeo local. 

= Separación de los soportes laterales para evitar la inestabilidad del perfil 
(más adelante se detalla este tema). 
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El momento plástico (Mp) de un perfil de acero es el momento que es capaz 
de plastificar todas y cada una de las fibras de la sección (en perfiles simétricos 
respecto al eje x-x la sección tiene la mitad de la sección en compresión y la otro 
parte en tracción, de modo que la sección se encuentre en equilibrio). Este valor 
de momento se encuentra relacionado con el módulo plástico (Z) que a su vez es 
la relación del área del perfil (A) multiplicado por distancia desde el centro de 
gravedad de cada área al eje neutro (d) como se presenta en la figura 13. 


| Zxx = A1 * d1 + A2 = d2 
a, 


eje neutro 


| Mp = Fy * Ly_x 
d, 





Figura 13 Esfuerzos internos en una sección sometida al momento plástico. 


Nota: Depende del criterio del diseñador elegir las secciones que cumplen 
con un momento máximo nominal relacionado con el momento elástico o plástico; 
de esta decisión dependerá el costo de la estructura y el nivel de seguridad, 
generalmente se intenta aprovechar al máximo las secciones y por ello se diseña 
para que el momento nominal sea Mn = Ø * Mp, donde Ø = 0.9; sin embargo, se 
debe asegurar las condiciones necesarias para llegar a la plastificación de las 
secciones. 


El programa Capacidad_Elementos_Acero_V2 determina la capacidad a 
carga axial, corte y momento de los elementos, posteriormente se habla con más 
detenimiento sobre esta función, por ahora se debe decir que la longitud no 
arriostrada de la viga debe ser menor que L, (para evitar pandeo local) y menor 
que L, para evitar pandeo lateral torsional, debe ser menor que las dos. 


3.2 CAPACIDAD AXIAL DE ELEMENTOS DE ACERO 


Al aplicar una carga axial de tracción en una sección de acero su 
capacidad nominal está dada por la ecuación 1. 


P, = F, * Az (1) 


Donde E, representa el esfuerzo de fluencia del material y A, el área de la 
sección transversal. No existe reducción de capacidad a tracción, no obstante, se 
limita la longitud del elemento hasta 300 veces el radio de curvatura mínimo de la 
sección, esto se lo hace para controlar flechas excesivas y ondulación o vibración 
lateral, además se procura que la rigidez sea adecuada con estos elementos. 
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Al aplicar cargas de compresión en un perfil de acero su capacidad se 
puede afectar debido a la inestabilidad o pandeo de la sección, estas fuerzas 
causan deflexiones laterales fuera de su plano y es por ello los códigos de diseño 
sísmico limitan la capacidad a compresión en función de un esfuerzo crítico 
(ecuación 2), que a su vez depende la relación de esbeltez 2 de la sección 
(ecuación 3). 


Pa = For * Ag (2) 


ja (3) 


Y 
Donde F„, es el esfuerzo critico; k, es un factor en función de las 
condiciones de borde o apoyo del elemento; L, es la longitud del elemento; r, es el 
radio de giro de la sección. Para el diseño con la metodología LRFD se debe 
reducir la capacidad nominal por factores @ que para compresión es 0.90 


La carga de pandeo de Euler es la carga axial máxima que se puede 
aplicar a una columna sin producir el pandeo de la misma. Esta carga es función 
del esfuerzo critico de Euler que se calcula con la ecuación 4. 


1? Es 


KL? (4) 
e 


FCYTEULER = 





Mientras mayor es la relación de esbeltez, hay mayor probabilidad de que 
la columna falle por pandeo y la sección tendrá menor capacidad. Por lo que la 
capacidad del material debe ser modificada por un factor, para determinar la 
capacidad real a la cual estará trabajando. Las ecuaciones que dispone la normas 
se realizaron en base a ensayos, en donde se evaluó experimentalmente la 
influencia de los esfuerzos residuales y las imperfecciones en los perfiles 
estructurales, generando una banda de posibles resultados como se muestra en la 
figura 14. 

OA 


Esbeltez media Esbeltez grande 
>< 





Banda de 
resultados 






Punto de 
inflexion 


¥ A 
Figura 14 Esfuerzo versus relación de esbeltez de un elemento sometido a 
compresión. 
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Existen varios autores que han propuesto curvas de esfuerzo — esbeltez, 
las cuales permiten determinar la capacidad real del material. 


La ecuación de Euler presenta valores superiores al límite de fluencia del 
material, es por ello que se corrige para calcular con ella el esfuerzo critico como 
se indica en la ecuación 5. 


FCTEULER < Fy > Fcr = FCTEULER 


(5) 


Fcreuer 2 Fy > Fcr = Fy 


La expresión de Rankine Gordon se presenta en la ecuación 6. 


Fcr x F 
Ferz EULER y (6) 
FCTeuLER + Fy 


La ecuación de Rankine es la ecuación 7. 
Fy 
2 
Fy + S (7) 


E 
it m** ES 


ror = 


La especificacion AISC 360 propone la ecuacion 8. 


Fcreuer < 0.44 : Fy > Fer = 0.877 * FCrEyLER 


Fy 
FCfrgyrgr = 0.44: Fy > Fer = 0.658F“"EULER x Fy 


Aplicando las expresiones para el material A36 se tiene en la figura 15. 
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Figura 15 Esfuerzo versus relación de esbeltez de un elemento sometido a 
compresión. 
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Para los cálculos se utiliza el valor dado por una de las curvas intermedias, 
es decir por la curva de la AISC 360, gráficamente se muestran en la figura 16. 


For \ m E s 


y 
e y A y 


Fy a \ 
E ASS SU 
0.44 Fy Sy” 0.877 Fi 
Inelástico Elástico 
4.71. [E a 
Fy 


Figura 16 Esfuerzo versus relación de esbeltez de un elemento sometido a 
compresión. 











4. PREDISEÑO DE ELEMENTOS 


Para el prediseño de elementos se asume que las vigas de los primeros 3 
pisos serán iguales (estructura de figura 6), independiente de que los vanos sean 
5 y 7 metros, y el diseño se realiza para la viga de mayor longitud, esto se lo hace 
por consideraciones prácticas constructivas y evitando que existan un desbalance 
con vigas en un mismo nudo. 


Inicialmente se recomienda determinar la altura del perfil, se sugiere tomar 
una altura mínima (h) en función de la longitud de la viga (L), como se recomienda 
en la ecuación 9. 


L 
=- (9) 
Para una longitud de 7 metros, la altura mínima de viga recomendada es 
0.40 metros; para definir las dimensiones de la sección se emplea la combinación 
de cargas más conservadora para cargas gravitacionales (NEC SE CG 2015), esto 
es: 


h 


C1 = 1.2: CM + 1.6- CV (10) 


La carga muerta CM en el piso inferior es 500 kg/m? y la carga viva CV es 
de 250 kg/m?. El ancho colaborante es 6 metros, por lo tanto, la carga 
uniformemente distribuida en la viga con la combinación de cargas se define por: 


q = (1.2: CM +1.6: CV) * Leoiaborante (1 1) 
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El valor de la carga “q” en la viga es de 6.0 T/m, con esta cantidad se 
determina los momentos en los extremos y en el centro de luz de la viga con las 
expresiones de la figura 17. Se obtiene un momento en el centro de 12.25 (T : m) y 
en los extremos de 24.5 (T : m). 





Figura 17 Diagramas de momento en una viga empotrada. 


Las vigas principales son parte del sistema de pórtico resistente a 
momento, es decir que además de soportar las cargas verticales, también deberá 
soportar cargas laterales, por lo tanto, el momento último sobre la sección debido 
a cargas verticales debe ser mayorado para considerar el efecto del sismo 
(Guerra, 2015); el valor del momento por sismo, para este caso se puede 
considerar como: 


Musismo = f. Mucarga gravitacional (1 2) 


La tabla 2 muestra los valores del factor f, para los distintos niveles de un 
edificio: 


Tabla 2 Factores de amplificación de efectos por sismo. 





Ya que se trata del prediseño para todos los pisos, se escoge el valor de 
f=1.2. Entonces el momento último de diseño, en su condición más 
conservadora es Musismo = 29.4 (T : m). 


Considerando el diseño de ultima resistencia donde se tiene que la 
capacidad del perfil (fp - Mn) debe ser al menos igual a la solicitación (Mu), para el 
prediseño se tiene: 

Ø: Mn ~ Mu 


Ø : Mp = Mu 


Mes (13) 


Ø 


Conocido el momento último, que se estima actuará en las vigas de Musismo = 
29.4 (T : m), se determina para este momento, el módulo plástico requerido: 
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Mp = Fy * Ly_x (14) 
Mu 
LXXreq = 0- Fy (15) 


Siendo el tipo de acero estructural A36, se tiene que Fy = 2530 [kg/cm’], 
ademas para el diseño a flexión el coeficiente de reducción es @ = 0.9, entonces: 





(29.4 Toma 1000 kg P 100 m) 


1T 1m 
0.9 « 2530-4, 
cm 


(16) 


LXXreq = 


(17) 
ZXxX req = 1291.17 [cm?*] 


En los catalogos comerciales se puede observar que la seccion IPE 450, 
cuyas dimensiones se aprecian en la figura 18, cumple con un modulo de sección 
elástico de 1500 cm’, es decir que su módulo plástico sera superior a este valor, 
por ello se considera esta sección para el análisis inicial. 


Datos de la sección 
IPE 450 





hs | 40o | mm | 
s= | 94 | m | 








Figura 18 Dimensiones de la sección IPE 450. 


Los autores consideran que este artículo debe ser una referencia para los 
estudiantes y personas que están iniciando en el calculo de estructuras, es por ello 
que se detallará el cálculo del módulo de sección elástico y plástico. 


Módulo de sección elástico: 


Para calcular el módulo elástico en una sección tipo “I” con patín superior 
e inferior de las mismas dimensiones se aplica la ecuación 18, la distancia ycg se 
mide desde el eje neutro a la fibra más alejada y es igual a h/2, además se 
determina la inercia del perfil con el siguiente razonamiento: 


Wy-x = Sy-x = (18) 





Nuevas Funciones de CEINCI-LAB para análisis y diseño... 17 










WS “ 
| 


m . 
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(b) 


Figura 19 a) Cálculo de la inercia de un perfil |. b) Resultados del cálculo. 
Módulo de sección plástico: 


En secciones simétricas se divide la sección como se indica en la figura 
20a y se determinar el módulo plástico de la sección con la ecuación 19, donde A; 
es el área de cada sección y d; es la distancia desde la mitad de la sección “l”, es 
decir la base de la “T”, hasta el centro de gravedad de cada sección (cada 
rectángulo). En la figura 20b se presenta los resultados del cálculo. 


Zxx = 2 » (A1 * d1 + A2 * d2) (19) 





(b) 


Figura 20 a) Cálculo del módulo plástico de sección. b) Resultados del cálculo. 
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Nota: Existe una diferencia entre el módulo presentado en el catálogo y el 
calculado debido a que las vigas IPE mantiene un radio de curvatura en la 
conexión del patín al alma, la variación es del 5% entre los valores mencionados. 


En la ecuación 20 se determina el factor de seguridad al utilizar esta 
sección, es decir la relación entre el módulo plástico de sección calculado de la 
viga IPE 450 versus el valor requerido. 


= LXXyiga _ 1623.9 Se (20) 
ZXXreg 1291.17 l 


Las secciones IPE generalmente cumplen con una compacidad para 
catalogar a las secciones como altamente dúctiles, esto se comprobará más 
adelante, pero de manera inicial se decide emplear estas dimensiones de vigas de 
los pisos 1 a 3 y para las vigas de la cubierta se empleará una sección IPE 400. 


Para prediseñar una columna se debe calcular la carga axial en función 
del área colaborante de la columna más cargada en planta baja, en este caso la 
columna del eje B de la figura 6. La expresión 21 permite realizar la sumatoria de 
la carga de cada piso. 


num pisos 


Pu = » (1.2: CMi + 1.6: CVi) [are i 


i=1 


Lyiga izquierda + Lyiga erect) (21 ) 
2 


El valor resultante es Pu = 128.9 (T). Amplificando este valor por f = 1.2 
debido a efectos sísmicos, como se indicó anteriormente para vigas, se tiene un 
valor de carga axial de PUsismo = 154.66 (T). 


Acorde a Guerra (2015) la relación de esbeltez para una columna de 
alrededor de 3 [m] de altura está entre 40 y 60. Para un elemento de acero A-36, 


con una relación de esbeltez de 60, el esfuerzo crítico será Fcr = 2093 H esto 
se indica en la tabla 3. 


Tabla 3 Esfuerzo critico calcula con diferentes metodologías. 
Datos del material 


2530.00 [Kg /enY] KL/r crítico = 
2.04E+06 [Kg /cm’] 


Fe EULER Euler corr Rankine Gor Rankine Fer AISC 360 
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Si se calcula la relación del esfuerzo determinado mediante las ecuaciones 
del AISC 360 y el Fy = 2530 [kg/cm2] del acero A36, se tiene una disminución 
del 20% de la capacidad del material. 


A continuación, se determina el área requerida de la columna como se 
indica en las siguientes ecuaciones: 








Fuerza (22) 
Esfuerzo = — 
Area 
, Fuerza PU Seine (23) 
BT CA requerida = Esfuerzo ~ @+ Fer 
154.66[T] - 000 ks 


j E 1T 
ATe€Qrequerida = 0.9 - 2093 [kg/cm?] 


AreQyequerida = 82.1 [cm?] 


Debido a cargas axiales se requiere una sección de columnas 
relativamente pequeña comparada con las vigas, no obstante, en columnas el 
diseño muchas veces no se encuentra controlado por fuerzas sino por los 
desplazamientos, específicamente para mantener derivas relativas de piso 
inferiores a las marcadas en los códigos de diseño, acorde con la NEC — 15 se 
tiene un límite máximo de 2%, valor que implica daño como se lo menciona en el 
artículo de Aguiar (2017 a). Incluso en muchos casos las dimensiones de las 
columnas se definen debido al diseño por capacidad, es decir que la capacidad de 
las columnas debe ser mayor a la de vigas para asegurar la disipación de energía 
como daño en estos elementos y esto se cumple cuando la relación indicada en la 
ecuación 24 es verdadera. 





ae > 1.0 (24) 


Donde Y Mpc es la sumatoria de los momentos plásticos nominales de las 
columnas que llegan a la junta; Y Mpv es la sumatoria de los momentos plásticos 
nominales de las vigas que llegan a la junta. Un diagrama de cuerpo libre de una 
junta o nudo para una columna interior de un pórtico a momento se muestra en la 
figura 21. 
































06 
Placas de Continuidad \ Mee Articulación plástica 
bi ! S | ] Mor 
| ~ | | | | 
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Figura 21 Nudo de una columna interior de un pórtico. NEC-SE-AC 2015. 
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Mediante un diseño en base a la metodología LRFD, es decir en base a 
factores de carga y resistencia, se define la sumatoria de momentos plásticos 
nominales de las columnas que llegan a la junta como se muestra en la ecuación 
25 y para vigas se presenta en la ecuación 26. 


X Mpc = Y Zot: (Fy 77) (25) 


> Mpv = ace - Ry - Fy : Zviga + Muyiga) (26) 


Donde Zcol y Zviga son los módulos plásticos de las columnas y vigas (si 
se considera una conexión RBS entonces Zviga debe calcularse en el centro de 
la RBS); Fy es el esfuerzo de fluencia del material que compone cada elemento; 
Puc es la carga axial requerida para la columna usando las combinaciones de 
carga para columnas (se debe incluir la sobre resistencia, más adelante se explica 
con detalle); Ag es el área bruta de la sección; Ry, es el factor de esfuerzo de 
fluencia probable (para A36 este valor es 1.3 acorde a la NEC-SE-AC 2015) y 
Muyigg €s el momento adicional basado en las combinaciones de carga 
(producido por la fuerza de corte localizado en la zona donde probablemente se 
tendrá la articulación plástica multiplicado por la distancia hasta la junta o nudo). 


Un índice para asegurar el criterio mencionado anteriormente, que se 
conoce como columna fuerte y viga débil, es la relación de la sumatoria de los 
módulos plásticos de sección de columnas respecto a la sumatorias de estos 
módulos de sección de vigas que debe ser superior a 1.15 como se indica en 
Uang et al (2001). A continuación, se presenta la deducción de esta relación 
considerando los parámetros particulares para este proyecto. 


Las columnas del eje B en el primer piso y segundo piso tendrán la misma 
sección, además se asume que la relación de la demanda versus la capacidad a 


oy . P 
compresión en la que se encontraran es del 30%, es decir que Jo" 0.3, 





Puc 


entonces I 0.3- Fy; por lo tanto, la ecuación 25 queda como se indica a 


Ym B Sl (r | 
pc = CO Ag 
X Mpc = 2:Zcol: (Fy — 0.3: Fy) 


X Mpc = 2-0.7-Zcol- Fy (27) 


continuacion. 


En las vigas si se considera una conexión RBS el módulo plástico de la 
sección disminuye debido a la reducción de los patines, por efectos de estar en un 
prediseño se considera una reducción de 30%, además por facilidad se puede 
asumir que Muyigg es el 15% del momento plástico de la sección (recordar que la 
viga izquierda y derecha tiene la misma sección). De la ecuación 26 se puede 
determinar que: 


> Mpv = X - Ry + Fy + Zviga_RBS + Muyiga) 
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> Mpv = 2: (1.1: 1.3 : Fy: 0.7  Zviga + 0.15: (1.1: 1.3 + Fy + 0.7: Zviga)) 
` Mpv = 2: 1.15 ; Fy : Zviga (28) 
Al reemplazar las ecuaciones 27 y 28 en 24 se tiene: 


Mpc 
RET, 
2, Mpv 





2 0.7 : Zcol : Fy x 
2: 1.15- Fy : Zviga — 





Zcol LGA 
Zviga (29) 


Para el caso particular la columna debe tener un valor de modulo plastico 
de sección de 1.64 veces el valor de Zviga = 1623.9 cm? esto implica un Zcol = 
2663.2 cm?. Además, se debe mencionar que en columnas se requiere de 
espesores de patin y alma superiores a los que generalmente se tiene en las IPE, 
es por ello que se recomiendan secciones tipo HEB. En este caso se determina 
una columna HEB 400 con un Zcol = 3125.4cm? y un área de 191.52 cm, 
además para las columnas de los ejes A y D se definen secciones HEB 360. 


Se debe mencionar que hasta este punto las secciones de los elementos 


presentadas en la figura 22 son en base a un prediseño y se requiere de cálculos 
que permitan corroborar las secciones de los elementos o refutar las mismas. 


IPE400 IPE400 IPE400 Story4 
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IPE450 IPE450 IPE450 
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Figura 22 Secciones del pórtico con base en un prediseño. 
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Las siguientes propiedades de material corresponden a las del acero 
ASTM A36, Fy = 36 ksi, Fu = 58 ksi y E = 29000 Ksi. Las dimensiones de las 
columnas y vigas se definen en la Tabla 4, donde la identificación del color es 
consistente con la figura 22. 


Tabla 4 Descripción de las secciones de los elementos estructurales. 
Color Elemento Dimensiones (mm) 





d x b x tf x tw 
| 
tf soon. 400X300x24x13.5 
Columnas 
tw q HEB360 360x300x22.5x12.5 
| beeen ABOX190x14.6x9.4 
TP Dees 400x180x13.5x8.6 





5. MODELACIÓN CON CEINCI-LAB 


Para comprobar las dimensiones de columnas y vigas se debe realizar el 
modelo del pórtico y en base al análisis sísmico determinar los desplazamientos, 
derivas y cortante basal, verificando que se cumpla con la deriva inelástica 
máxima de la NEC-15 para estructuras de acero. Además, se debe definir las 
combinaciones de carga y determinar las fuerzas internas en los elementos, para 
pórticos planos esto implica Axial, Corte y Momento Flector. Posteriormente se 
debe comparar la demanda máxima con la capacidad de los elementos para 
verificar que sean adecuados, además se comprueba que la relación de la 
capacidad de las vigas versus la capacidad de las columnas que convergen a un 
nudo sea inferior a 1, asegurando un comportamiento sismorresistente adecuado. 


Se debe mencionar que las unidades en CEINCI-LAB deben ser 
consistentes durante todo el análisis, por ejemplo, en ese caso se emplea 
Toneladas fuerzas y metros. 


Inicialmente se requiere de datos los vectores sv, sp que corresponden a 
las separaciones de cada vano y alturas de piso; con las funciones 
geometria_nudo_viga; glinea_portico2; gn_portico2, con ello se calcula la 
cantidad de vanos, pisos, columnas, vigas, nudos y restricciones, además se 
definen las coordenadas X, Y de la estructura y los vectores NI, NJ, que contienen 
a los nudos iniciales y finales de los elementos. 


En el programa se considera cargas uniformemente distribuidas en cada 
piso (Carga por unidad de longitud), esto se asigna en los vectores CM y CV que 
corresponden a los estados de carga muerta y viva. Para el análisis sísmico se 
considera que la carga reactiva del sismo se compone por la carga muerta, en el 
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caso general, y casos particulares se considera un porcentaje de la carga viva, 
para este caso se asume un 25% de la carga viva adicional a la carga muerta total. 


Para el análisis sísmico de pórticos planos, considerando únicamente la 
componente horizontal del sismo, se puede definir un modelo simplificado de un 
grado de libertad por piso como se indica en Pilatasig et al. (2020). Para la 
definición de este modelo se calculan los grados de libertad con la función 
cg_sismo2. El código de CEINCI-LAB hasta este punto del análisis se muestra en 
la figura 23. 


Código del programa principal, parte 1: 


%% Geometría del pórtico 

sv =[5.00;7.00;5.00]; %Ingresar un vector con la longitud de vanos en 
(metros) 

Long_Piso=5*2+7.0; % Longitud (m) 

Ancho_Col=6; % Ancho colaborante 

sp =[3.6;2.88;2.88;2.88]; %lngresar un vector con la altura de pisos en 
(metros) 


[nv,np,nudt,nudcol,nudvg,nod,nr|=geometria_nudo_viga(sv,sp); 
[X, Y]=glinea_portico2(nv,np,sv,sp,nod,nr); 
[NI,NJ]=gn_portico2(nr, nv, nudt, nudcol, nudvg); 


%% Cargas aplicadas a las vigas de cada piso 
CM=[0.5;0.5;0.5;0.35]*Ancho_Col; % Carga muerta T/m en cada piso 
CV=[0.25;0.25;0.25;0.1]*Ancho_Col; % Carga viva T/m en cada piso 
Porc_viva=0.25; % Fracción de carga viva para la carga reactiva del sismo 
(W=D+Porc*L) 

for i=1:np 

qp(i,1)=-(CM(i,1)+(Porc_viva*CV(i,1))); % Carga distribuida en vigas en cada 
piso T/m 

end 

[CG,ngl]=cg_sismo2(nod,nr, Y); %Calcula los grados de libertad análisis 
sísmico 


% % Dibujo de la Estructura 2D 
% dibujo(X,Y,NI,NJ) 
dibujogdl(X, Y,NI,NJ,CG) 
dibujoNudElem(X, Y,NI,NJ,CG) 





Figura 23 Codigo de CEINCI-LAB con datos de geometria y cargas. 


Con las rutinas de CEINCI-LAB dibujo, dibujogdl, dibujoNudElem se 
puede graficar el esquema estructural; los grados de libertad en cada nudo y la 
numeracion de nudos y elementos, como se puede observar en la figura 24a y 
24b. Los grados de libertad presentados son para el analisis sismico. 
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Esquema Estructural 





(b) 
Figura 24 a) Grados de libertad para análisis sísmico; b) Numeración de nudos y 
elementos. 


En la monografía de Aguiar et al. (2020) se indica el procedimiento para 
asignar diferentes secciones a los elementos del pórtico, esto es en base a la 
matriz Sec_VG_COL que contiene en cada fila un tipo de sección, además en la 
primera columna de cada fila se indica el tipo de sección seguido de las 
dimensiones del elemento; las dimensiones son peralte, d; espesor del alma, tw; 
ancho del patín superior, bs e inferior bi; espesor de los patines superior tfs e 
inferior tfi. 


Para designar cada sección a un elemento se debe tener en cuenta la 
numeración de los mismos, que se indica en la figura 24b, la matriz 
Elem_Sec_VG_COL que tiene en la primera columna el tipo de sección y seguido 
a esta columna se debe escribir el número de elemento a la cual se asignara la 
misma. El código de CEINCI-LAB se expone en la figura 25. 
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Código del programa principal, parte 2: 


SECCIONES 
% Propiedades del material 
Fy=25310.505; % Esfuerzo de fluencia, T/m2 
Es=20389019.16;% Módulo de elasticidad, T/m2 
Fu=40778.036; % Esfuerzo último, T/m2 
Uso Conexion = 1; 
Ry = 1.3; % Factor de Fluencia Probable (Ry) depende del acero para A36 
Ry=1.3 


%lUso Conexion = # % Conexión Precalificada RBS 1, Conexión 
Precalificada diferente que RBS 2, No Precalificada 3. 
Materiales = [Fy, Fu, Es, Uso_Conexion]; 


% Sec_VG_COL = Elementos | =[tipo#, d,tw,bs,bi,tfs,tfi]... En milimetros!! 


sec_VG_COL=[ 1 400, 13.5, 300, 300, 24.0, 24.0; % Tipo 1 HEB 400 
2 360, 12.5, 300, 300, 22.5, 22.5; % Tipo 2 HEB 360 
3450, 9.4, 190, 190, 14.6, 14.6; % Tipo 3 IPE 450 
4 400, 8.6, 180, 180, 13.5, 13.5]; % Tipo 4 IPE 400 


% Elem_Sec_VG_COL = ["tipo# respecto de Sec_VG_COL", Identificador de 

los elementos] 

Elem_Sec VG_COL=|1 2, 6, 10, 14, O; 
1 3; 7,11,15, 
21, 5, 9,13, 0, O; 
24, 8,12,16, 0, 0; 
3 17, 18,19, 20, 21,22; 
3 23, 24, 25, 26, 27, 28; 
3 29, 30, 31, 32, 33, 34; 
4 35, 36, 37, 38, 39, 40]; 


, 0; 
0 





Figura 25 Código de CEINCI-LAB con datos de materiales y secciones. 


En este punto se tiene la geometría, materiales y secciones del pórtico en 
análisis, en primera instancia se procede con el análisis sísmico asumiendo que se 
encuentra emplazado en Quito, en suelo tipo D. Aplicando la NEC-SE-DS, se 
obtiene los siguientes parámetros que intervienen en el cálculo del porcentaje de 
la carga sísmica que constituirá el corte basal acorde al “Método Estático 


Equivalente” (tabla 5); se tiene una columna que indica la referencia de la NEC- 
SE-DS 2015. 
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Tabla 5 Parámetros que intervienen en el cálculo del Corte Basal. 
CARGAS DE SISMO 


Zona Sísmica Tabla 1. Sec.3.1.1 Pág. 27 
Factor de aceleración en la zona Tabla 1. Sec.3.1.1 Pág. 27 
Tipo del perfil del suelo Tabla 2. Sec.3.1.1 Pág. 30 
factor de sitio Fa Tabla 3. Sec.3.2.2 Pág. 31 
factor de sitio Fd . Tabla 4. Sec.3.2.2 Pág. 31 
factor comportamiento inelás. suelo i Tabla 5. Sec.3.2.2 Pág. 32 
Factor asociado al periodo de retorno Sec 3.3.1 Pág. 34 

Relación de amplificación espectral : Sec.3.3.1 Pág. 34 


Aceleración en T=To l S, =NZFa Sec.3.3.1 Pág. 34 


pan F 
Periodo Límite en T=To Ty = 0.10F, 2 Sec.3.3.1 Pág. 35 


Periodo Límite en T=Tc . r.=0.5SF..4 Sec.3.3.1 Pag. 34 
c= 
Periodo Limite en T=TL Sec.3.3.1 Pag. 34 


Coeficiente Ct Sec.6.3.3 Pág. 62 
Coeficiente para Calculo de Periodo l = 2.4.F Sec.6.3.3 Pag. 62 


Altura total del elemento . Planos 
Período teórico método 1 Sec.6.3.3 Pág. 62 





Período teórico método 1 mayorado Sec.6.3.3 Pág. 63 


Donde V es el cortante basal; Sa(T) es la aceleración espectral que se 
determina en función del periodo fundamental de vibración de la estructura T como 
se indica en la ecuación 30; R es el factor de modificación de respuesta sísmica; 
Op y De es el factor de irregularidad en planta y elevación; W es la carga reactiva 
al sismo. En la tabla 6 se presenta los valores de las variables mencionadas y el 
calculo del coeficiente sísmico, el cortante basal es 19.84% de la carga reactiva 
total. Se debe mencionar que el valor de R se considera como 6, en la normativa 
se permite un valor máximo de 8. 


Sa =nZFa para O<T<Tc 


Tex” (30) 
Sa = nZFa (=) para T >Tc 


Donde y es un factor que depende de la región; Z es un factor de la 
aceleración de la zona; Fa, es el coeficiente de amplificación de suelo en la zona 
de periodo corto; Tc, representa el periodo límite de la meseta en el espectro de 
aceleraciones de la NEC 15. 


En este punto es importante destacar que los espectros de diseño, en 
general, para un período igual a cero parten del PGA y en un período T, alcanzan 
la meseta del espectro definida por la primera ecuación de (30). En forma 
conservadora no se ha considerado esta rama del espectro es por eso que se 
tiene la primera rama entre 0 <T <Tc. 


Los períodos que definen el punto inicial y final de la meseta del espectro 
que son Tọ y T,, están en función de los factores de sitio de aceleración F,, de 
desplazamiento F, y de suelo F, estipulados en NEC-15. Ahora bien, para la 
ciudad de Quito, es conveniente trabajar con los factores de sitio que obtuvo ERN 
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(2012) en el estudio de microzonificación sísmica de la ciudad. Para facilitar su uso 
se han elaborado dos versiones de MIZOSIQ que es el programa que, en forma 
amigable, reporta estos factores para cualquier lugar de la ciudad de Quito. 
Sangurima (2013); Astudillo (2014); Aguiar (2017 b). 


Tabla 6 Coeficiente sísmico. 


Pórtico de acero 


Parámetro Variable Valor Unidad Obs. Referencia 
Aceleración en T=To Sao T11<TC Sec.3.3.1 Pág. 34 
Factor de irregularidad en planta Op 1 - Tabla 13. Sec. 5.2.3 Pag. 50 
Factor de irregularidad en 













elevacion ee 1 - Tabla 14. Sec. 5.2.3 Pag. 51 
Factor de reducción de R 6 

respuesta - Tabla 16. Sec. 6.3.4 Pág. 65 
Factor de importancia | 1 - Tabla 6. Secc. 4.1 Pág. 39 
Coeficiente sísmico Coef. 0.1984 - - Sec. 6.3.2 Pág. 61 
AO On k 1.0575 - T>05 Sec. 6.3.5 Pag. 67 


fuerzas laterales 


En la NEC-SE-DS 2015 se indica el procedimiento para determinar las 
fuerzas laterales en cada piso, para ello se aplica la ecuación 31 que se muestra a 
continuación: 


Wjx*hj* 


"J = Sn Wis hit 


(31) 


Donde Wj, Wi es la fracción del peso en cada nivel; hj, hi es la altura de 
los niveles en j e i; el valor de “k” se relaciona con el periodo de la estructura y 
varía entre 1 y 2, este valor refleja la forma en la que se distribuyen las fuerzas 
laterales en cada nivel de la estructura; es decir que la distribución es lineal 
(triangular) para valores de T < 0.5 segundos como se muestra en la figura 26, 
mientras que para estructuras de mayor altura (con periodos mayores a 0.5 
segundos) este valor se calcula como k = 0.75 + 0.5 T y sik = 2 la distribución 
es parabólica (el valor máximo de k es 2); V es el cortante basal. 





En gs War n a Min 
Te Š 
Fi Wi 1 rm 
> T $ 
k T * 
FAN Wi | o m 
hi 
© Je | J $ 
Fy Gal | Wi 1 ¿m 
k ee ¥ oe. 
$$ $$ 
Vi 


Figura 26 Fuerzas laterales en cada piso. 


Otra metodologia para calcular el cortante basal es el Método Modal 
Espectral, que se basa en discretizar a la estructura en sus modos de vibración 
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como se muestra en la figura 27 (formas de vibrar), a cada modo se asocia una 
frecuencia y una forma modal. Mediante el espectro se evalúa la respuesta de 
cada modo de vibración de la estructura y posteriormente se obtiene la máxima 
respuesta global de la estructura en términos de desplazamientos, derivas y 
cortantes en cada piso, para ello se combinan los efectos de cada modo. 



































e l re | 
MODO | MODO 2 MODO 3 


Figura 27 Ejemplo de los modos de vibración de un pórtico plano. 


La combinación modal que se emplea en CEINCI-LAB es denominada 
CQC, también conocida como Combinación Cuadrática Completa, que considera 
el amortiguamiento modal como se muestra en la ecuación 32. 





La combinación CQC se usa para encontrar los desplazamientos laterales 
y para las fuerzas laterales, destacando que en el último caso el criterio se aplica 
directamente en las fuerzas laterales en cada modo. En Aguiar (2012) se 
presentan dos formas de encontrar las fuerzas, la primera es a partir de los 
cortantes y la segunda como se hace en este artículo. 


Para realizar el análisis sísmico con el Método Estático Equivalente y 
Modal  Espectral se desarrolla una nueva función denominada 
Analisis Sismico Porticos; esta función internamente realiza el calculo de la 
matriz de rigidez del pórtico, de la matriz de masas y calcula los periodos, las 
frecuencias y modos de vibración. Se determina el periodo con la ecuación de la 
NEC 15 (Tygc = 0.534 seg se indica en la tabla 5), se verifica que el periodo de 
vibración calculado con la masa y rigidez sea menor a 1.3 veces el periodo de la 
NEC 15, caso contrario continua el cálculo con el periodo de la NEC mayorado en 
1.3; este es un control importante para evitar que las fuerzas sísmicas sean 
reducidas debido a la flexibilidad de la estructura. Se recomienda que el periodo 
de la estructura no supere este límite de periodo. 


La función de CEINCI-LAB también permite calcular el cortante basal, la 
distribución de las fuerzas laterales en la altura del pórtico, desplazamientos y 
derivas inelasticas de cada piso (ecuación 33). 
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Am = 1.0-Ae-R (33) 


Despl. Rel. 
Ae = | (1) 


(34) 


Donde Am es la deriva inelastica para cada piso, este valor debe ser 
menor a 0.02 para pórticos de acero (NEC 2015); R es el factor de modificación de 
la respuesta sísmica; Ae es la deriva elástica para cada piso; Despl. Rel.) es el 
desplazamiento de cada piso obtenido con la aplicación de fuerzas laterales; h es 
la altura de cada piso. 


El código de CEINCI-LAB se describe en la figura 28. 


Código del programa principal, parte 3: 


%% Cálculos Sísmicos 

% Espectro de diseño NEC 15 

is=4; % Suelo D 

iz=5; % Zona sísmica 5 = 0.4g 

lp=2; % Región Sierra 

Importancia = 1.0; % Factor de importancia 

R = 6; % Factor de reducción de la fuerza sísmica 

Cd = 5.5; % ASCE 7 -16 SMF 

omega=3; % Sobre resistencia 

Dibujar = 1; % 1 para dibujar el espectro // 2 no dibujar 

fip=1 ;fie=1.0; % Factores de planta y elevación 

zeda=0.05; % Amortiguamiento 

Tipo_Estr =1; 

% Tipo_Estr =1 Acero sin arriostramientos // Tipo_Estr =2 Acero con 
arriostramiento 

% Tipo_Estr =3 Hormigón armado sin arriostramientos // Tipo_Estr =4 
Hormigón armado con arriostramiento 

Espectro=[is iz ip Importancia R fip fie zeda Tipo_Estr Dibujar]; % Datos para 
el espectro NEC 15 


[M_Masa, M_Rigidez Sis, T, T_Nec, Cort_Sis, Fuerzas Sis, Derivas, 
Desp_Inelast]=Analisis Sismico Porticos(sv, sp, qp, Sec_VG_COL, 
Elem_Sec_VG_COL, Materiales, Espectro); 


disp('Los periodos de vibración en segundos con CEINCI-LAB son:’) 
T% seg 


disp('El periodo de vibración fundamental en segundos NEC 15 es:') 
T_Nec % seg 





Figura 28 Código para análisis sísmico con CEINCI-LAB. 
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La matriz de rigidez se muestra a continuación, se puede observar que la 
dimensión de esta matriz es de 4x4 debido a que se considera un modelo 
simplificado de 1 grado de libertad por piso, entonces la matriz de rigidez de la 
estructura con los grados de libertad indicados en la figura 24a se condensa en los 
grados de libertad principales (los horizontales): 


22480 -16252 4414.3 -543.96 
-16252 25923 -16443 3458.2 
4414.3 -16443 224/76 -9899.8 
-543.96 3458.2 -9899.8 6909.1 


La matriz de masa también tiene dimensión 4x4 debido a que el modelo 
concentra la masa del piso en un punto, es decir una masa en cada nivel. 


5.85 0 0 0 
M= 0 9.85 0 0 
0 0 9.85 0 

0 0 0 3.90 


El periodo fundamental de vibración calculado con la matriz de masa y 
rigidez es de 0.615 segundos, que es inferior al limite de 1.3 veces el periodo de la 
NEC 15. En la figura 29 se puede visualizar el espectro elástico de diseño de color 
azul y de color verde el espectro reducido por el factor R, además la línea vertical 
de color rojo representa el periodo fundamental de vibración y las otras líneas de 
color gris son los otros periodos de vibración. 


nl Espectro de Diseño NEC 15, Amortiguamiento 5% 





Espectro Elástico 
Espectro Inelástico R = 6.00 
—---T fundamental = 0.62 seg 





Aceleracion (g) 








0 0.5 1 1.5 2 2.5 3 3.5 4 
Periodo (seg) 


Figura 29 Espectro de diseno. 


En la figura 30 se muestra los resultados de desplazamientos, derivas y 
cortante en cada piso del portico, determinados con ambas metodologias, la linea 
de color rojo representa el Método Estatico Equivalente y la de color azul el 
Método Modal Espectral. En la figura 30 se tiene las respuestas mediante los 
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métodos de análisis mencionados, generalmente las respuestas del análisis Modal 
Espectral son inferiores a las respuestas del método Estático Equivalente. El 
usuario debe escoger con que fuerzas laterales continua el análisis, no obstante, 
se debe cuidar que la respuesta del análisis Modal Espectral sea superior al 80% 
del método Estático Equivalente si el usuario decide emplearlo en los siguientes 
análisis. En este caso el cortante basal con el método Estático Equivalente es 41.7 
T y con el método Modal Espectral es 36.6 T; esto representa 87.8 % entre un 
método y otro. Además, la deriva inelastica máxima que se presenta en el 
segundo piso es 1.34 % y 1.16 % con cada método, en ambos casos es menor al 
2% que se tiene como límite en la NEC 15 para pórticos de acero. Cuando se 
considera una conexión RBS se debe amplificar las derivas para considerar el 
cambio de área e inercia en las secciones de vigas, se recomienda amplificar 
estos valores por 10% para ser conservadores y verificar que sean menores al 
límite máximo de las normas, en este caso se cumple. 


DESPLAZAMIENTO DERIVA POR PISO CORTE POR PISO 
T T T T 

















Altura de pisos 
Altura de pisos 


Altura de pisos 
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Figura 30 Respuestas del análisis sísmico. 


A continuación, se debe calcular el índice de estabilidad de la estructura 
(9), esto se relaciona con el efecto P-Delta, es decir con efectos adicionales que 
pueden causar un incremento en las fuerzas internas, momentos y derivas, esto se 
debe considerar para la evaluación de la estabilidad estructural ante posible volteo 
y se determina con la ecuación 35. 





0; = —— (35) 


Donde Pi es la carga axial debido a los estados de cargas muerta y viva, 
es posible considerar la reducción de cargas vivas como se establece en las 
normativas, sin embargo, para ser conservadores en CEINCI-LAB no se considera 
reducción. A, es el desplazamiento relativo elástico de cada piso, si se considera 
una conexión RBS debido a la reducción de los patines en las vigas se puede 
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aumentar un 10% estos desplazamientos; V; es el cortante en cada nivel; h; es la 
altura de cada nivel. En la figura 31 se presenta el código de CEINCI-LAB, la 
función Indice _Estabilidad_Sismo_V2 utiliza para el cálculo de la estabilidad las 
fuerzas y desplazamientos obtenidas con el método Estático Equivalente, esto se 
define en la variable Metodo que es igual a 2 en este caso, si se requiere calcular 
con el método modal espectral esta variable debería definirse como 1. 


Código del 


% Análisis de estabilidad 
Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Estático Equivalente Nec 15 


[Teta, Teta2, Teta_Admisible]=Indice_Estabilidad_Sismo_V2(CM, CV, 
Long_ Piso, sp, R, Metodo, Cort_Sis, Desp_Inelast, Uso Conexion) %NEC 15 





Figura 31 Código para análisis de estabilidad con CEINCI-LAB. 


Si en el análisis se incluye los efectos P-Delta desde un principio, 
entonces se puede calcular el índice de estabilidad con la ecuación 36. 





(36) 


Cuando el índice de estabilidad es menor a 0.10 entonces se clasifican los 
pisos como no susceptibles al volteo; si este índice es superior a 0.3 entonces la 
estructura es potencialmente inestable y debe rigidizarse; cuando se encuentra 
entre 0.1 y 0.3 se debe encontrar un factor fp_, y multiplicar todas las fuerzas 
laterales por este factor. 


1 


fp-a = 1-08 (37) 


Siendo @ el mayor valor de los 0,. A continuación, se muestran los 
resultados del índice de estabilidad de cada piso de la estructura y se comprueba 
que es menor a 0.10 en todos los niveles, por lo tanto, se puede continuar el 
análisis sin consideraciones adicionales. 


0.0176 
0.0180 
0.0118 
0.0059 


Para corroborar de manera parcial que las secciones de columnas y vigas 
son adecuadas se procede a calcular la relación de los módulos plásticos de 
columnas respecto a vigas, que llegan a un mismo nudo, asegurando que sea 
mayor a 1.64 como se demostró en la ecuación 29. La rutina que calcula esta 
relación se denomina relacion_sum_Z columna_viga y 
dibujo relacion_sum_Z columna_viga permite presentar gráficamente esta 
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relación (figura 32), es importante mencionar que los códigos sísmicos indican que 
en el ultimo piso no se requiere cumplir con esta relación de la capacidad de las 
columnas superior a la capacidad de las vigas. 


Código del 


[L,seno,coseno]=longitud(X,Y,NI,NJ); %long. de (Velem, seno y coseno 

[Seccion,Elem_Tipo Prop2,Z ELEM,th_ELEM]=Diferentes Vg COL(nudt,Sec_VG 
_COL,Elem_Sec_VG_COL,Uso_Conexion,Fy,Es,L); 

%Seccion=[ELEMENTO, Area, Inercia, NUMEROS DE ELEMENTOS A CREAR, 
INCREMENTO EN NUMERO DE ELEMENTO] 


[ELEM]=gelem_portico(Seccion); 


Capacidad_sum_Z columna_viga=relacion_sum_Z columna_viga(coseno,nr 
snudcol,X,Y,NI,NJ,Z ELEM,Elem_ Tipo _Prop2,Sec_VG_COL); 


dibujo relacion_sum_Z columna_viga(X,Y,NI,NJ,Capacidad_sum_Z column 
a_viga) 





(a) 


Verificación de Capacidad *suma Z Columna/Viga 


olVg = 1.3 ColNVg = 2.1 





ol g = 1.9 Col/Vg = 3.2 


6 = 
ol/Vg = 3.2 ol/Vg = 1.9 ol/Vg = 1.9 ColWVg = 3.2 
4+t 


2F ColWVg=3.2 olg = 1.9 olg = 1.9 Colg = 3.2 














(b) 
Figura 32 a) Código de CEINCI-LAB; b) Relaciones de módulos plásticos de 
columnas para vigas de forma gráfica. 


En esta etapa del análisis se comprueba que las secciones de columnas y 
vigas asignadas al pórtico permiten cumplir con los requisitos “viga débil y columna 
fuerte”, estabilidad global y derivas, por lo tanto, se puede continuar con el análisis 
estático, es decir aplicar los estados de carga gravitacionales y laterales (sismo) 
para determinar las fuerzas internas en los elementos. 


Para el análisis estático se cambian los grados de libertad y se considera 
un modelo con 3 grados de libertad por nudo (Desplazamiento en horizontal, 
vertical y Giro), esto debido a que es un pórtico plano. En la figura 33a se muestra 
el código para calcular los grados de libertad con la rutina cg, es importante notar 
que se requiere definir las restricciones en la base de cada columna, para este 
caso se considera un empotramiento. En la figura 33b se muestra la numeración 
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de grados de libertad con la función dibujogdl. Con este modelo se calcula la 
matriz de rigidez KTT. 


Código del incipal, parte 6: 


%% Análisis Estático 
aux=length(sv); 
for i=1:aux+1 
RES(i,:)=[i 1 1 1]; % Restricciones de apoyo 
end 
[CG,ngl]=cg(nod,nr, RES); %Calcula los grados de libertad 


dibujogdl(X,Y,NI,NJ,CG) 
[VC]=vc(NI,NJ,CG); %Vector de colocación 
%% Matriz de Rigidez 


[K]=krigidez_acero(ngl, ELEM,L,seno,coseno,VC,Es); 
KTT=K; %< MATRIZ "K" DE TODA LA ESTRUCTURA (T/m) 


(a) 


Esquema Estructural 








(b) 
Figura 33 a) Código de CEINCI-LAB; b) Grados de libertad para el análisis 
estático. 


Se deja al lector ejecutar el programa hasta este punto y visualizar la 
matriz de rigidez, se puede mencionar que la matriz tiene dimensiones de 84x84 
debido a que es una matriz cuadrada en función de los grados de libertad, que 
para este caso es 84. Aguiar (2020) 


A continuación, se aplica las cargas sísmicas en la estructura, como se 
indica en la figura 34, recordando que un caso de carga es cuando se analiza con 
fuerzas laterales de izquierda a derecha y otro caso es aplicando las fuerzas en 
sentido opuesto, esto es importante en estructuras que no son simétricas. 
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Figura 34 Cargas laterales aplicadas a las estructuras. 


Código del incipal, parte 7: 


%% Cargas Sísmicas 

Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Estático Equivalente Nec 15 

njc=np; % Numero de juntas Cargadas 

% Sentido de izquierda a derecha 

F = [[nr+1:nr*2-1:nod]' Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % F = 
[Nudo cargado, FH, FV, Momento] 

datos=0; 

nmc=0; % Numero de miembros cargados 

Fm=0; %Elem carg, carga, codigo, elem a gener, incr numero elemento 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 

% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 

% Desplazamientos y Giros 

q_Sismo=KTTIQ_Sismo; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 

[Fuerzas_Elementos Sismo_P]=fuerzas_acero_arm1_V2(ngl,ELEM,L,seno,cosen 


0,VC,Es,q_Sismo,Q2_Sismo,CG,nr,NI,NJ,nmiembro,RES); 


% Sentido de derecha a izquierda 

F = [[nr+2*nr-1:nr*2-1:nod]' -Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % F = 
[Nudo cargado, FH, FV, Momento] 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 

% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 

% Desplazamientos y Giros 

q_Sismo=KTTIQ_Sismo; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
[Fuerzas_Elementos_ Sismo _N)]=fuerzas_acero_arm1_V2(ngl,ELEM,L,seno,cosen 
o,VC,Es,q_Sismo,Q2 Sismo,CG,nr,NI,NJ,nmiembro,RES)); 


Fuerzas Elementos Sismof[1,1)=Fuerzas Elementos Sismo P; 
Fuerzas Elementos Sismof[2,1)=Fuerzas Elementos Sismo_N; 





Figura 35 Código de CEINCI-LAB. 
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En la figura 35 se muestra el código de CEINCI-LAB que permite realizar 
el análisis estático con fuerzas laterales, el procedimiento consiste en definir las 
fuerzas en los nudos, esto se lo hace en el vector F. Se calcula el vector de cargas 
generalizadas, es decir las cargas en los nudos de las estructuras que son 
correspondientes a los grados de libertad indicados en la figura 33 b. Se resuelve 
la operación matricial planteada en la ecuación 38, donde Q y K son el vector de 
cargas generalizadas y la matriz de rigidez de la estructura (Aguiar 2020). 
Entonces se calculan los deslazamientos y giros del vector q (se debe considerar 
los grados de libertad de la figura 33 b). 


Q=K-q (38) 
Codigo del incipal, parte 8: 


%% Cargas gravitacionales 
njc=0; % Numero de juntas Cargadas 
F = 0; % F = [Nudo cargado, FH, FV, Momento] 
datos=0; 
nmc=nudvg; % Numero de miembros cargados 
ul=1; 
for i=1:np 

for j=1:nv*2 
Fm_CM(ui,:)=[nudcol+ui CM(i,1) 1 0 0]; %Elem carg, carga, código, elem a 
gener, incr numero elemento 
Fm_CV(ui,:)=[nudcol+ui CV(i,1) 1 0 0]; %Elem carg, carga, código, elem a 
gener, incr numero elemento 
ui=ui+1; 

end 
end 
[Q_CM,Q2_CM]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm_CM,datos, 
nr, NI,NJ); 
[Q_CV,Q2_CV]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm_CV,datos,nr 
,NI,NJ); 
% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 
% Desplazamientos y Giros 
q_CM=KTTQ_CM; 
q_CV=KTTQ_ CV; 
% Fuerzas y momentos finales en los elementos 
nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
[Fuerzas Elementos CM]=fuerzas_acero_arm1_V2(ngl, ELEM,L,seno,coseno,VC 
,Es,q_CM,Q2_CM,CG,nr,NI,NJ,nmiembro,RES)); 


[Fuerzas Elementos CV]=fuerzas_acero_arm1_V2(ngl, ELEM,L,seno,coseno,VC, 
Es,q_CV,Q2 _CV,CG,nr,NI,NJ,nmiembro,RES)); 





Figura 36 Código de CEINCI-LAB. 


Con el vector de coordenadas generalizadas q y la matriz de rigidez de 
cada elemento k se calcula las fuerzas en cada elemento, esto se lo hace en 


Nuevas Funciones de CEINCI-LAB para análisis y diseño... 37 





coordenada locales; cada elemento tiene 6 grados de libertad, 3 en el nudo inicial 
y 3 en el nudo final. La matriz de celdas Fuerzas Elementos Sismo contiene en 
su posición {1,1} las fuerzas en los elementos debido a las fuerzas laterales de 
izquierda a derecha y en la posición {2,1} las fuerzas en los elementos debido a 
fuerzas laterales en sentido opuesto. 


En la figura 36 se muestra el código que permite calcular las fuerzas en los 
elementos debido a caga muerta y viva por separado, los resultados se almacenan 
en Fuerzas Elementos CM y Fuerzas _ Elementos CV respectivamente. El 
procedimiento es similar al indicado anteriormente, con la particularidad de que las 
fuerzas en los nudos no existen, entonces la variable F se define por cero y como 
se tiene cargas uniformemente distribuidas en las vigas, estas se definen en las 
matrices Fm_CM y Fm_CV para carga muerta y viva. 


Los estados de carga se pueden superponer para considerar las 
combinaciones de cargas mencionadas en los códigos de diseño, gracias a este 
principio los efectos de cargas en los elementos se suman directamente entre los 
diferentes estados de carga. 


Para el diseño de los elementos acorde al método LRFD, se debe 
incrementar las solicitaciones en base a unas combinaciones de cargas; en la 
NEC-SE-CG 2015 se tienen las siguientes combinaciones de carga: 


Combinación 1 1.4D 

Combinación 2 1.2D + 1.6L + 0.5max|[Lr; S; R] 
Combinación 3 1.2D + 1.6max[L; S; R] + max[Lr; 0.5W] 
Combinación 4 1.2D + 1.0W + L + max![Lr; S; R] 
Combinación 5 1.2D + 1.0E + L + 0.2S 
Combinación 6 0.9D + 1.0W 
Combinación 7 0.9D + 1.0E 


Donde D corresponde a la carga muerta (carga permanente + peso 
propio); E carga de sismo (componente horizontal y vertical); L sobrecarga viva; Lr 
es la sobrecarga viva en cubierta; S carga de granizo; R es la carga de lluvia; W 
carga de viento. 


En las combinaciones de carga donde se incluye los efectos del sismo es 
posible hacer una reducción del efecto de la carga viva en función de la 
probabilidad de ocupación, no obstante, no se considera esta reducción para ser 
conservadores en el análisis. También se debe considerar que las combinaciones 
de carga con sismo implican la componente vertical, pero en estos cálculos no se 
incluye este efecto. Con estas consideraciones se calcula las combinaciones de 
carga aplicables en este caso particular, con los efectos de carga muerta, viva y 
sismo (con fuerzas laterales en ambos sentidos); en la matriz Combinaciones se 
tiene las fuerzas en los elementos con las combinaciones de estados de carga y 
en la matriz Envolvente se determina los máximos efectos de estas 
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combinaciones en cada elemento (máximos y mínimos). El código de CEINCI-LAB 
se presenta en la figura 37, se indica que la función se denomina 
Combinacion_Cargas. 


Código del incipal, parte 9: 


%% Combinaciones de casos de cargas 


[Combinaciones, Envolvente] = Combinacion Cargas (Fuerzas Elementos CM, 
Fuerzas Elementos CV, Fuerzas _ Elementos Sismo); 





Figura 37 Código de CEINCI-LAB. 


En las combinaciones de carga que involucran los efectos del sismo debe 
considerarse la sobre resistencia (factor 2), esto se aplica para comprobar el 
diseño de elementos que deben estar sujetos a deformaciones plásticas y que 
puedan comprometer la estabilidad de la estructura, por ejemplo, en la verificación 
de la capacidad a compresión de columnas. 


En la sección 4.2 de la NEC-SE-AC 2015, se menciona la utilización del 
factor de sobre resistencia en las combinaciones de carga con las que se diseña 
ciertos elementos y conexiones del sistema de cargas laterales, no obstante, no se 
brinda mayor detalle de como se debe considerar este factor y es por ello que se 
acogen las directrices de las normativas internacionales como el código ASCE 7 
16 (Minimum Design Loads and Associated Criteria for Buildings and Other 
Structures). 


En la sección 2.3.6 del ASCE 7-16, se indica que cuando una estructura 
este sujeta a efectos sísmicos entonces se debe considerar también el efecto de 
combinaciones de cargas adicionales con el efecto de sobre resistencia y en su 
sección 12.4.3.1 se muestra que este factor amplifica el efecto de la componente 
horizontal del sismo, como se indica en la ecuación 35. 


Emn = 2° Qe (39) 


Donde Emn es la componente horizontal del sismo con el efecto de sobre 
resistencia; R es el factor de sobre resistencia (en el ASCE 7-16 se exponen 
valores referenciales en función del sistema estructural); Qg es el efecto de sismo 
horizontal (cortante basal). 


Asumiendo que se desprecia el efecto de sismo vertical entonces las 
combinaciones de carga 5 y 7 cambian como se indica a continuación. 


Combinación 5b 1.2D + 1.0E mn + L + 0.25 


Combinación 7b 0.9D + 1.0E mh 


Para considerar estos nuevos estados de carga se hace un análisis similar 
al caso anterior, con la particularidad de que se debe mayorar las fuerzas laterales 
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debido al sismo en su componente horizontal, el código que genera esto se 
muestra en la figura 38. Las variables que consideran esta modificación se 
presentan con la distinción de la palabra omega y de color rojo. Este factor fue 
definido anteriormente con las propiedades del sistema estructural y toma un valor 
de 3. 


Código del incipal, parte 10: 


%% Cargas Sísmicas CON FACTOR DE SOBRE RESISTENCIA 

Metodo = 2; % 1 Modal Espectral Nec 15 // 2 Lineal Estático Nec 15 

njc=np; % Numero de juntas Cargadas 

% Sentido de izquierda a derecha 

F = [[nr+1:nr*2-1:nod] omega*Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % 
F = [Nudo cargado, FH, FV, Momento] 

datos=0; 

nmc=0; % Numero de miembros cargados 

Fm=0; %Elem carg, carga, codigo, elem a gener, incr numero elemento 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 

% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 

% Desplazamientos y Giros 

q_Sismo=KTTIQ_Sismo; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
[Fuerzas_Elementos Sismo _P_ omega]=fuerzas_acero_arm1_V2(ngl,ELEM,L,sen 
o,coseno,VC,Es,q_Sismo,Q2_ Sismo,CG,nr,NI,NJ,nmiembro,RES); 


% Sentido de derecha a izquierda 

F = [[nr+1:nr*2-1:nod]' -omega*Fuerzas_Sis(:,Metodo) zeros(np,1) zeros(np,1)]; % 
F = [Nudo cargado, FH, FV, Momento] 
[Q_Sismo,Q2_Sismo]=cargas_arm(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos, 
nr,NI,NJ); 

% [Q,Q2]=cargas(njc,nmc,ngl,L,seno,coseno,CG,VC,F,Fm,datos); 

% Desplazamientos y Giros 

q_Sismo=KTTIQ_ Sismo; 

% Fuerzas y momentos finales en los elementos 

nmiembro= 1; % Número del Primer elemento a partir del cual se calcula la 
contribución a la matriz de rigidez 
[Fuerzas_Elementos _Sismo_N omega]=fuerzas_acero_arm1_V2(ngl,ELEM,L,sen 
o,coseno,VC,Es,q_Sismo,Q2_Sismo,CG,nr,NI,NJ,nmiembro,RES); 


Fuerzas Elementos _Sismo_omegaf[1,1)=Fuerzas Elementos Sismo _P omega; 
Fuerzas _Elementos_Sismo_omega[2,1)=Fuerzas Elementos Sismo _N_ omega; 


%% Combinaciones de cargas CON FACTOR DE SOBRE RESISTENCIA 
[Combinaciones_omega, Envolvente_omega] = 
Combinacion_Cargas(Fuerzas_Elementos_CM, Fuerzas_Elementos CV, 
Fuerzas_Elementos_Sismo_omega); 





Figura 38 Código de CEINCI-LAB. 
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En la sección 6.2 de la NEC-SE-AC 2015 se manifiesta que columnas y 
vigas de pórticos especiales a momento deben tener secciones clasificadas como 
sismicamente compactas, esto se puede interpretar como secciones de alta 
ductilidad como lo indica el AISC 341-16 en su capitulo D. Este requerimiento es 
importante debido a que asegura que las secciones no sean susceptibles de fallas 
anticipadas por inestabilidad local o pandeo local. 


Tabla 7 Limites de ancho versus espesor para elementos sometidos a 
compresión. Fragmento de la tabla D1.1 del AISC 341-16. 


Límites de 
A 


Referencia 
Ana Amd 
Elemento de Alta Elemento de Moderada 
ductilidad ductilidad 


Para Ca < 0.114 Para Ca < 0.114 


E E 
29 1 — 1.04 .96 1 — 3.04C€ 
57 0 104c,) | 3 Em a) 


Para C, > 0.114 Para C, > 0.114 


E E 
0.88 Dea | Cee) 
Rr g) y: y a 


2 1.57 


Py = Ry: Fy: Ag Py = Ry: Fy: Ag 





En la tabla D1.1. del AISC 341-16 se muestran los límites para clasificar a 
las secciones según la relación de aspecto de patines y alma denominada A; si 
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este factor es menor a Aa la sección es altamente dúctil y si se encuentra entre 
Ana Y Ama Se clasifica como moderadamente dúctil; si el valor de A es superior a 
este límite entonces la sección es no compacta. 

Para secciones tipo “lI” los patines son secciones no atiesados o 
rigidizados y se debe probar con el caso 1 de la tabla D1.1 y para el alma, como 
es un elemento atiesado entonces se debe verificar en función del caso 7 de la 
misma tabla en el AISC 341. En la tabla 7 se presentan estos casos de 
verificación. 


Código del incipal, parte 11: 


%% Compacidad de elementos considerando SMF * AISC -341-16 
[Compacidad, Lambda, relacion_aspecto]=Compacidad_ Elementos(nudcol, 


nudt,Es,Fy,Ry,ELEM,Elem_Tipo_Prop2,Sec_VG_COL,Envolvente_omega); 


dibujo Compacidad(X,Y,NI,NJ, Compacidad) 
% HD=Altamente ductil // MD = Moderadamente Ductil // ND = No Ductil 


(a) 


Esquema Estructural: Compacidad 











(b) 
Figura 39 a) Codigo de CEINCI-LAB; b) Clasificacion de las secciones en el 
portico. 


El Cálculo se realiza con la función Compacidad_Elementos indicada en 
la figura 39a, como resultados se tiene la matriz compacidad que almacena la 
clasificación del patín, del alma y de la sección de cada elemento con códigos 1 
para HD, 2 para MD y O para ND, en este caso HD significa que la sección es 
altamente dúctil y se tuviera la nomenclatura MD es moderadamente dúctil, 
finalmente para secciones con poca ductilidad se establece la nomenclatura ND. 
La matriz Lamba contiene los límites para HD o MD de patín y alma de cada 
elemento. Finalmente, en la matriz denominada relacion_aspecto se almacena la 
relación de ancho para espesor de patín y alma de cada elemento. Además, se 
presenta gráficamente los resultados con la función dibujo_Compacidad como se 
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muestra en la figura 39b, los códigos HD, MD y ND también definen un color al 
elemento, como azul, verde y rojo respectivamente. Para este caso todas las 
secciones son altamente dúctiles. 


A continuación, se detalla el cálculo para la clasificación de la sección 
correspondiente a la columna del primer piso el eje B, el perfil es un HEB400 
cuyas dimensiones son: peralte total d = 400 mm; ancho del patín bt = 300 mm; 
espesor del alma tw = 13.5 mm; espesor del patín tf = 24 mm. 


Para empezar, se debe realizar la verificación de cada elemento por 
separado, es decir del patín y del alma de la sección. 


e Para el patín: 


ON 20389019.16 
oe = Fy 1325310505 
20389019.16 
Ama = 0.4: 
= Fy 1.3 - 25310.505 _ 


“| 








Como la sección “I” es simétrica y para ser conservadores el valor de b se 
puede tomar como la mitad del ancho del patín. 


b _bt/2 _ 300/2 
=" = 6.25 
tf f 24 


Como 1 = 6.25 es menor que Apa = 7.96 entonces el patín es altamente 
ductil (HD). 


e Para el alma: 


Se debe calcular la relación Ca, en ella interviene la carga axial última Pu, 
se debe indicar que este valor se obtiene de la envolvente con las combinaciones 
de carga considerando el factor de sobre resistencia en el cálculo y el factor @ se 
considera 0.9 para ser conservadores. En la figura 40 se puede apreciar que entre 
mayor sea el factor Ca entonces los límites para clasificar a una sección como 
altamente dúctil o moderadamente dúctil son más estrictos, es decir que son más 
pequeños. 


La columna es el elemento número 2 en CEINCI-LAB (figura 24b), 
entonces se puede hallar la carga axial última con el factor de sobre resistencia 
como el valor numérico correspondiente al mayor valor entre 
Envolvente_omega{1,1}(2,1) y Envolvente_omega{1,2}(2,1) sin considerar el 
signo, es decir valor absoluto. Para este caso particular se tiene Pu = 139.557. 


Py = Ry-Fy- Ag =1.3-25310.505 - 0.019152 = 630.17 T 
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e a Pu _ 18955 e 
7 @.:Py 0.9-630.17 ` 


Límites para clasificar un elemento como HD o MD 
120 





0 


0.2 0.4 0.6 0.8 1 
Ca 


——)\ hd ——Amd ----- Ca=0.114 





Figura 40 Límites para clasificar una sección HD o MD. 


Como Ca es mayor que 0.114 entonces se debe calcular los limites Ana y 
Ama con las siguientes ecuaciones: 


Ang = 0.88 : - (2.68 — C,) = 0.88 : oa - (2.68 — C,) 
M = Fy 1.3 -25310.505 a 
Ang = 0.88 alee ded (2.68 — 0.246) = 53.32 
eo 1.3-25310.505 ~ a 
eS ee 20389019.16 ES 
Fy Ca) = 1.3 -25310.505 A 

Ama = 1.29 A (2.12 — 0.246) = 60.18 
ee 1.3-25310.505 `“ OS 


Se debe comprobar que estos valores sean superiores al siguiente límite 


LE 20389019.16 4. og 
| = Fy 1.3 -25310.505 


Se determina la relación de aspecto como 





F 


Ama = 1.29 





f 





h d—2-tf 400—2-24 
LS Ss 507 
tw tw 13.5 


Como A = 26.07 es menor que Ara = 53.32 entonces el alma es altamente 
dúctil (HD). 
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e Enconclusion: 


Como los patines y el alma se clasifican como altamente dúctiles entonces 
la sección se clasifica como altamente dúctil (HD). 


Para el análisis de los otros elementos la metodología es similar, en vigas 
se podría asumir que C, = 0 debido a que las solicitaciones axiales en estos 
elementos son muy bajas. 


Los elementos de un pórtico plano están sometidos a cargas axiales, de 
corte y a momentos flectores, la función que permite calcular la capacidad nominal 
de los elementos incluyendo la reducción por los factores @ se denomina 
Capacidad Elementos Acero (figura 41). Esta función calcula la capacidad 
nominal reducida por axial, corte y momento flector en los vectores phi_Pn; 
phi_Vn; phi_Mn; además calcula el vector Longitudes que contiene los valores 
de Lp y Lb que deben ser corroboradas por el usuario para colocar 
arriostramientos laterales a los elementos y con ello evitar el pandeo lateral 
torsional y pandeo local, esto es importante para que el elemento pueda 
desarrollar su máxima capacidad a flexión. El vector Longitudes tiene una 
particularidad y es que solo calcula los valores de Lp y Lb para cada tipo de 
secciones, por ejemplo, en este pórtico se tiene 4 tipos de secciones de 
elementos, que se definieron en la variable Sec VG_COL (figura 25). 
Factor K Long Efectiva calcula los factores k para determinar la longitud 
efectiva de los elementos sometidos a compresión. 


Código del incipal, parte 12: 


%% Capacidad de elementos de acero 
[Ki]=Factor_K Long _eEfectiva(coseno,nr,nudcol,X,Y,NI,NJ,L,ELEM); 


[phi_Pn,phi_Vn,phi_Mn,Longitudes]=Capacidad_Elementos_ Acero _V3(Es,Fy,Ry,Z 
_ELEM,Seccion,L,Ki,Elem_Tipo_Prop2,Sec_VG_COL, coseno); 
% Arriostrar a los elementos con una longitud menor al mínimo entre (Lp, Lb) 





Figura 41 Código de CEINCI-LAB. 


La metodología para calcular la capacidad axial de un elemento se 
presentó en la sección 3.2 de este artículo. En este apartado se muestra el 
procedimiento para determinar el factor K, el mismo depende de las condiciones 
de borde del elemento, como se muestra en la figura 42. 


En pórticos de acero una forma para conocer el factor K para cada 
columna es empleando los nomogramas propuestos por Jackson y Moreland como 
se muestra en la figura 43. Los factores G se determinan con la ecuación 40, y los 
subíndices A y B se refieren a base y cabeza de columnas. 
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La geometría deformada 
se muestra en línea 
discontinua 


Valor teórico de K 


Valor recomendado 
para condiciones 
aproximadas 





Figura 42 Factores K para diferentes condiciones de apoyo. 
360 — 10. Alacero. 





0.65 








0.80 


WHT 











THVT 





2.0 


y 2.0 


Tomado de: AISC 





Figura 43 Nomograma para determinar los Factores K en marcos rigidos que 


permiten el desplazamiento lateral. Tomado de: AISC 360 — 10. 
Alacero. 


CT), 


Gamm 
E), 


Donde I es la inercia de la sección en sentido x-x; L es la longitud del 


(40) 


elemento; los subíndices c y b se refieren a columnas y vigas respectivamente. 
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Para la programación de estos monogramas se pueden emplear la 
ecuación 41, que tiene como solución aproximada la denominada regla francesa 
indicada en la expresión 42. 


Gar Gor (y) -36 (E) _ 
6- (Ca + GR) tan() 





Gat Ga +75 


Si el extremo de la columna tiene un apoyo articulado entonces se puede 
tomar un valor de G = 10, si tiene un empotramiento G = 1 en ese borde. 


Se detalla el calculo para la columna en planta baja del eje B, en la figura 
44 se muestra los valores de G, y Gg en los extremos de la columna, además se 
presentan los valores de inercia y longitud en cada elemento que rodea al nudo 
superior (B). 


Ic, = 0.000558710784 m* 
Lc, = 2.88 m 


© 
IPE450 G, =317 IPE450 


Iv, = 0.0003214043879744 m' Iv, = 0.0003214043879744 m“ 
Lv, =5.0m Lv, = 7.0m 
HEB400 
Ic, = 0.000558710784 m‘ 
Lc, = 3.6m 


0 


Figura 44 Detalle de valores involucrados en el cálculo de K para la columna de 
planta baja en el eje B. 





Como en vigas y columnas se tiene el mismo material entonces el valor de 
Gb se determina de la siguiente manera: 


Ic, y 10 
_ Le, beg | 
Ce = To, Iv T T 
Lv, Lv, 


El factor K para esa columna es: 
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E 1.6- Ga’ Gp +4- (Ga +G) +7.55 _ 16:1:317+4 (1 +3.17) +75 g 
7 Ga + Gg +7.5 7 1+3.7+7.5 o 


Para la columna se tiene un factor K = 1.58 


También se podría utilizar el nomograma presentado en la figura 43, 
donde se debe considerar como un punto en la recta numérica izquierda el valor 
de G, y un punto en la recta derecha el valor de Gg, luego al unir estos puntos 
mediante una recta entonces la intersección de la misma con la recta numérica del 
centro indicará el valor del factor K. Esto se ilustra en la figura 45, obteniendo un 
factor K de 1.58. 








Figura 45 Determinación del factor K mediante el nomograma para marcos rígidos 
que permiten el desplazamiento lateral. 


El valor de esfuerzo crítico de Euler se determina como: 


1? Es 1? Es T? -20389019.16 [T 
e E E O (=) 


DD 
Tex K-L 1.58 : 3.6 
a [0.000558710784 
A 0.0192 
Calculando con base en la especificación AISC 360, presentada en la 


ecuación 8. Como FCcrtyrer-xx > 0.44: Fy entonces la capacidad axial nominal a 
compresión se calcula como 





Fy 25310.505 T 
Fez = 0.658FTEULER-ax + Fy = 0,65818099105 - 25310.505 = 2.3872 - 10* (7) 


Se verifica que Fcr,, es Menor que Fy. 
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Pnyx = FCT ex, * A = 2.3872 - 104 - 0.0192 = 458.34 (T) 


La capacidad nominal a compresión se ve afectada por = 0.9, por lo 
tanto: 
Ø- Pn, = 412.51 (T) 


Es importante mencionar que este calculo se realiza en el plano de 
analisis, no obstante, la columna puede pandear fuera del plano; debido a que la 
columna es tipo “I” este segundo escenario (pandeo fuera del plano) tiene mayor 
probabilidad de pasar porque es el plano débil de la sección. Por lo tanto se 
calculará la capacidad axial en sentido transversal, para simplificar el análisis se 
asume que el factor K es igual al caso anterior, en un análisis espacial se debe 
determinar analiticamente este valor en función de las vigas que convergen al 
nudo en sentido perpendicular al plano XY. 


El valor de esfuerzo crítico de Euler se determina como: 


1? Es 1? Es T? -20389019.16 


T 
Pomos == = == = == = 3.5009 + 10% (=) 





Kiyo 
a K-L 1.58 - 3.6 


Ie [0.00010807 
A 0.0192 
Calculando con base en la especificación AISC 360, presentada en la 


ecuación 8. Como FCcrtyrer-xx = 0.44: Fy entonces la capacidad axial nominal a 
compresión se calcula como 


Fy 25310.505 T 
Fcr,, = 0.658F¢TEULER-zz - Fy = 0.6583.5009-104 - 25310.505 = 1.8702 : 10* (—) 
Se verifica que Fcr,, es menor que Fy. 
Pn,, = Fcr,, : A = 1.8702 : 104 - 0.0192 = 359.078 (T) 
La capacidad nominal a compresión se ve afectada por = 0.9, por lo 
tanto: 


Ó-Pn,, = 323.17 (T) 


La capacidad a compresión de la columna debe ser la menor entre Ø. 
Pn,, y Ø: Pn,,, para este caso la capacidad axial es: 


Ø- Pn = 323.17 (T) 


En vigas generalmente el efecto de la carga axial no es influyente, estos 
casos no se calcula el factor k, sino que es posible asumir como k = 1. 
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La capacidad nominal a corte (Vn) se determina con la ecuación 43. 
Vn =0.6:Fy: Aw: Cv, (43) 


Donde Aw es el área del elemento que resiste el corte, en perfiles tipo | se 
asume que Aw es el área del alma y se calcula con la ecuación 44; Cv, es un 
factor debido a la resistencia del corte y se determina con la ecuación 45. 





Aw =(d-—2:tf):tw (44) 
= Cuando kis <1.1: pa 
tw Fy 
Cv, = 1 
= Caso contrario 
A (45) 
e ET 


El valor de la variable kv depende de la existencia de rigidizadores 
transversales; si no se tiene estos elementos entonces kv = 5.34 que es el caso 
conservador con el que se calcula en CEINCI-LAB. Cuando existan rigidizadores 
transversales entonces kv se debe calcular con la ecuación 46. 


ky =5+ Ea 
(r) (46) 


k, = 5.34 cuando a > 3.0 


Donde a es la distancia libre entre rigidizadores. 


A continuación, se ejemplifica el cálculo de la capacidad a corte de las 


columnas HEB400. 
h _d-2 tf 400 2:24_ | 
tw tw 135 © 


i kok a 5.34 + 20389019.16 _ ee 
Fy ` 25310.505 a 


h Ky E : 
Como Jab s entonces se considera Cv, = 1 





Aw = (400 — 2:24) : 13.5 = 4752 mm? = 4752: 10"? (m?) 
Vn = 0.6: Fy: Aw: Cv, =0.6:25310.505 - 4752 - 1076 - 1 = 72.165 (T) 
La capacidad nominal a corte se ve afectada por pọ = 0.9, por lo tanto: 


Ø- Vn = 64.95 (T) 
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En la sección 3 del presente artículo se realiza una introducción al 
comportamiento de elementos sometidos a flexión, se indica que la capacidad de 
los elementos depende de la compacidad de la sección para evitar pandeo local, 
no obstante, también depende de la longitud del elemento para controlar el posible 
efecto de pandeo lateral, es decir de una inestabilidad del elemento fuera del 
plano de aplicación de cargas (ver figura 46). 


La figura 46 muestra la reducción de capacidad en función de la longitud 
no arriostrada de los elementos, es decir que una sección altamente dúctil con una 
longitud no arriostrada menor a Lp podrá desarrollar la máxima capacidad a 
flexión, pero si su longitud es mayor a Lp esta capacidad disminuirá. 

















PANDEO PANDEO PANDEO 
FLEXIONAI PORSIONAI FLEXOTORSIONAI 


Figura 46 Pandeo de un elemento en un plano distinto al de aplicación de cargas. 


A 





= 
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Resistencia Nominal de 
Flexión, Mn 
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Lp Lr 
Longitud No Arriostrada, Lb 


Figura 47 Capacidad a flexión en función de la longitud no arriostrada el elemento. 


Las longitudes Lp, Lr se pueden calcular con las ecuaciones 47 y 48. 


E 
Lp = 1.76: ry: 5 (47) 
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Md << (48) 
r JST 0.7 Fy 
ly:C 
rż = a — (49) 
Para secciones H con doble simetria y alas rectangulares entonces: 
ly +h? (50) 
E 


Donde ry es el radio de giro de la sección respecto al eje y-y, en este 
caso corresponde al eje débil de las secciones l; J, es la constante torsional de la 
sección; Sx es el módulo elástico de la sección respecto al eje x-x; Para secciones 
| con simetría c = 1. ho es la distancia entre los centroides de los patines. 


La capacidad a flexión de una sección | con doble simetría y compacta se 
puede determinar con la ecuación 51. 
Mp =Fy:Z, Cuando Ly < Ly 
Lb — Lp 
Mn =; Cb: [Mp — (Mp — 0.7: Fy: Sx): E] <Mp Cuando Lp<Lb<Lr (51) 
Fcr Sx < Mp Cuando Lb > Lr 


Donde Cb es el factor de modificación por pandeo lateral-torsional para 
diagramas de momento no uniformes y Fer es un esfuerzo crítico reducido para 
flexión, las expresiones matemáticas para calcular estos valores se encuentran en 
el capítulo F del AISC 360 -16. 


También se debe mencionar que para que una sección sea considerada 


altamente dúctil se debe tener un arriostramiento del elemento a una distancia no 
mayor a Lb, que se calcula con la ecuación 52. 


Lb = 0.095 -ry 





By Fy (52) 


En la función Capacidad Elementos Acero V3 de CEINCI-LAB la 
capacidad a flexión se calcula con la ecuación 53, pero para que esto se cumpla 
en la realidad se debe verificar que las secciones sean altamente dúctiles y que 
estén arriostradas a una distancia menor entre el mínimo de Lb y Lp. 


Mn = Mp = Fy : Z; (53) 


Anteriormente se detalló el cálculo del módulo plástico de sección para la 
columna HEB400, cuyo valor es Zcol = 3125.4 cm? considerando la ecuación 19, 
además se demostró que es una sección altamente dúctil. En este apartado se 
procede a calcular los límites de arriostramiento Lb y Lp. En primer lugar, se 
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requiere el radio de giro ry y la inercia Iy que se pueden calcular con las 
ecuaciones 54 y 55. 








ly 
ey 124 54 
ry A (54) 
a, fan B-72) o (DAY ge 
y=(5)-2 — OA (85) 


El área de la sección transversal es 0.0192 mY; la inercia Ix es 0.0005587 
m* y la inercia en el eje débil de la sección, respecto al cual el elemento es 
susceptible de pandeo lateral, ly es 0.00010807 m’. 


_ [ly _ [000010807 _ _... 
= 1A I opie O 
E 


A AA 
Ry-Fy ` 073 "7 3-25310505 A m) 


Lp = 1.76 is 1.76 - 0.075 pues 3.75 
ee ae y arene 25310505.) 


Las columnas con sección HEB400 requieren tener una longitud sin 
arriostramiento lateral menor a 3.75 (m), que es el mínimo valor entre Lb y Lp. En 
este caso las columnas de primer piso tienen una longitud de 3.6 metros, en ese 
nivel las columnas se encuentran arriostradas por las vigas y losas del sentido 
transversal, por lo tanto, es válido calcular la capacidad a flexión de esta sección 
con la ecuación 53. 


Lb = 0.095 -ry : 








Mn=Mp=Fy:Z = 25310.505 (<) 3125.4 =) ( ma ) 
cia as lea aa ee Nae A cr 


Mn = 79.106 (T : m) 
La capacidad nominal a flexion se ve afectada por = 0.9, por lo tanto: 
$: Mn = 71.2 (T:m) 
Para las columnas HEB360 los límites Lp y Lb son 3.81 y 4.49 m 
respectivamente, por lo tanto, al tener una longitud sin arriostrar menor a estos 


límites entonces es aplicable la ecuación 53 para hallar su capacidad a flexión. 


En las ecuaciones 47 y 48 se puede apreciar que, para un mismo tipo de 
acero, las longitudes Lp y Lb solo dependen del radio de giro ry. 


En la figura 48 se muestra la variación de estos valores para diferentes 
secciones IPE, representadas por su radio de giro ry. En esta figura se aprecia 
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que las vigas IPE450 requieren de arrostramientos laterales separados por una 
distancia máxima de 2 metros entre ellos, una forma de brindar el arriostramiento 
lateral en vigas del pórtico resistente a momento es mediante las vigas 
secundarias colocadas en estos intervalos. 


Si se cumple con una distancia no arriostrada menor a Lp para vigas, 
entonces su capacidad a flexión se determina con la ecuación 53, esto lo calcula 
automáticamente CEINCI-LAB, no obstante, si no se cumple con ello entonces el 
usuario debe calcular la capacidad de estos elementos de acuerdo a la ecuación 
51, considerando como Lb la longitud real no arriostrada. 


O 
ow 
ona. 
3 9.6 ty Y 
o O4 a 
on O,, Wa 
b OM Wg * @ 
= aa w £” 0 o 
oOo NN => o 
£ > oo N y o > e. 
O ON W e 
-l O Oru a. ~ a o 
15 oN~ wat © o 
1 waa 
egoe 
2 B @Lb (m) eLp (m) 
0.5 ? 


0 1 2 3 4 5 
Radio de giro, ry (cm) 


Figura 48 Longitud máxima sin arriostramiento para vigas IPE (acero A36). 


Para la comprobación de la demanda versus la capacidad de los 
elementos, en base a la metodología de Factores de Carga y Resistencia (LRFD 
por sus siglas en inglés, Load and Resistance Factor Design) se debe cumplir que: 


D- Ra >=) Vi Qi = Qu (56) 


Donde @ son factores de reducción de la capacidad nominal Rn; y; son 
factores que amplifican los efectos de las cargas; Q; corresponde a las cargas 
aplicadas; Q,,, es la carga última o efecto en los elementos con las solicitaciones 
amplificadas. 


Es importante destacar que la capacidad axial y a flexión calculadas son 
valores máximos que pueden resistir los elementos, no obstante, estos efectos se 
encuentran interrelacionados como se puede observar en las ecuaciones 57 y 58. 





Pu Pu 8 Mu 
<1 (57) 


P S09 as 2 o oe 
“2 hs PA ®-Pn 9 @-Mn 
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u P Pu ii Mu Ži (58) 
l > — + — 
Ø : Pn 2:0:Pn Ø- -Mn 





Para 


En la figura 49 se puede apreciar el diagrama de interacción entre el 
efecto axial y flexión para una columna W8x31, donde se puede entender como la 
capacidad a axial se reduce a medida que la se tiene solicitaciones a flexión. 






Curva Promedio Aproximada 
para Perfiles H 


Carga Axial Normalizada, P/P, 


0.0 0.2 0.4 0.6 0.8 1.0 
Momento de Flexión Normalizado, Mpc /M, 


Figura 49 Diagrama de interacción Axial-Flexión. Tomado de: AISC 360 — 10. 
Alacero. 


En CEINCI-LAB se han calculado las solicitaciones con las combinaciones 
de carga de la NEC-15 y también la capacidad de los elementos, por lo tanto, se 
puede calcular la relación entre la demanda versus la capacidad y verificar que 
estos factores sean menores a 1 para asegurar que las dimensiones de los 
elementos sean adecuadas para resistir las solicitaciones. Para la verificación de 
la capacidad a compresión en columnas mediante la ecuación 59, el usuario debe 
conocer que la fuerza axial última considera el efecto de la sobre resistencia para 
amplificar la componente horizontal del sismo. Para comprobar el resto de 
elementos o incluso las columnas ante otros efectos de solicitaciones, las 
combinaciones de carga no incluyen la sobre resistencia. 








Pu 

Ø : P, > Pu > a (59) 
Vu 
Mu 

> < 1 

Ø: M, > Mu 5. Mn ` (61) 


En la figura 50a se presenta el código para determinar la demanda versus 
la capacidad a compresión de las columnas, se destaca que las solicitaciones 
incluyen el factor de sobre resistencia (Q). En la figura 50b se presenta los 
resultados de forma gráfica y se observa que la relación D/C es menor a 1. 
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Código del incipal, parte 13: 


%% Verificacion de la capacidad axial de las columnas con sobre resistencia 

for i=1:nudcol 
D_C_Axial_omega(i,1)=max([abs(Envolvente_omega{1,1}(i,1)); 

abs(Envolvente_omega{1,2}(i,1))])/ phi_Pn(1,1); % Nudo inicial 


D_C_Axial_omega(i,2)=max([abs(Envolvente_omega{1,1}(i,4)); 
abs(Envolvente_omega{1,2}(i,4))])/ phi_Pn(1,1); % Nudo final 
End 


dibujo Demanda _Capacidad_Axial Columnas Sobre resistencia(X,Y,NI,NJ, 
D_C Axial_omega); 





(a) 


Verificación de Demanda/Capacidad Axial con sobre resistencia 


D/C = 0.044 D/C = 0.064 D/C = 0.064 
D/C = 0.14 D/C = 0.18 D/C = 0.18 
D/C = 0.24 D/C = 0.3 D/C = 0.3 





(b) 
Figura 50 a) Codigo de CEINCI-LAB; b) Relaciones demanda/capacidad de 
columnas incluyendo el efecto de sobre resistencia. 


En la figura 51 se presenta el código para generar las comprobaciones a 
flexo-axial, para pórticos de acero generalmente los efectos axiales son más 
críticos a compresión y por ello se ha decidido denominar efectos de flexo 
compresión. También se verifica la relación demanda versus capacidad a corte de 
los elementos. Las funciones denominada dibujo _*---” son las encargadas de 
representar las relaciones demanda versus capacidad como se indican en las 
figuras 52a y 52b, los elementos adoptan un color en función de la relación D/C, 
es decir que se mantiene de color cyan cuando D/C es menor a 0.5; mayor a 0.5 
cambia a color verde hasta que D/C sea 0.7; luego se convierte en amarillo hasta 
una relación D/C de 0.9 y magenta hasta un valor de 0.95, finalmente al pasar este 
limite los elementos adoptan el color rojo. 
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Código del incipal, parte 14: 


%% El resto de elementos se verifica sin sobre resistencia 
for i=1 :nudt 

D_C_Axial(i,1)=max([abs(Envolvente(1,1)(1,1)); aobs(Envolvente{1 ,2}(i,1))])/ 
phi_Pn(i,1); % Nudo inicial 

D_C Axial(i,2)=max([abs(Envolvente{1,1}(i,4)); aos(Envolvente(1,2)(1,4))])/ 
phi_Pn(1,1); % Nudo final 


D_C Corte(i,1)=max([abs(Envolvente({1,1}(i,2)); aos(Envolvente(1,2)(1,2))])/ 
phi_Vn(1,1); % Nudo inicial 

D GC Corte(i,2)=max([abs(Envolvente{1,1}(i,5)); abs(Envolvente/1,2)(1,5))])/ 
phi_Vn(1,1); % Nudo final 


D_C_Flexion(i,1)=max([abs(Envolvente{1,1}(i,3)); abs(Envolvente{1 ,2}(i,3))])/ 
phi_Mn(1,1); % Nudo inicial 
D_C_Flexion(i,2)=max([abs(Envolvente{1,1}(i,6)); abs(Envolvente{1 ,2}(i,6))])/ 


phi_Mn(i,1); % Nudo final 


a tinh cae a 1),D_C_Axial(i,2)); 
D_C_Vn_ví(i,1)=max(D_C_Corte(i,1),D_C_Corte(i,2)); 
D_C _Mn_vf(i,1)=max(D_C_Flexion(1,1)D_C_Flexion(i,2)); 


fD GC Pn ví(i,1)<0.2 
D GC Flexo Compr(i,1) =D_C Pn_vf(i,1)/2+D_C_Mn_ví(i,1); 
else 
D_C_Flexo_Compr(i,1) = D_C_Pn_vf(i,1)+8/9*D_C_Mn_vf(i,1); 
end 
end 


%% Dibujos de Demanda/Capacidad 
dibujo Demanda_Capacidad_Flexo Compresion(X,Y,NI,NJ,D_C_Flexo_Compr); 
dibujo Demanda_Capacidad_Cortante(X,Y,NI,NJ,D_C_Vn_vf); 





Figura 51 Código de CEINCI-LAB 


En las figuras 52a y 52b se puede observar que las relaciones de 
demanda vs capacidad de los elementos no sobrepasan a 1, es decir que las 
secciones impuestas en el prediseño cumplen con la resistencia ante las 
solicitaciones de forma satisfactoria. 


Como paso final se debe verificar que la relación de la capacidad de las 
vigas versus la capacidad de las columnas sea inferior a 1 (ecuación 62) de esta 
forma se tendrá que las columnas son más fuertes que las vigas y la disipación de 
energía se dará en las vigas. 


Y Mpv 


smpe £ (62) 
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Verificación de Demanda/Capacidad *Flexo-Compresión 





D/C = 0.57 D/C = 0.79 D/C = 0.79 D/C = 0.57 


D/C = 0.11 


D/C = 0.16 


D/C = 0.19 





D/C = 0.19 D/C = 0.21 D/C = 0.21 D/C = 0.19 


0.0 | 0.9 | 0.95 
(b) 


Figura 52 a) Demanda vs Capacidad a Flexo-Compresión b) Demanda vs 
Capacidad a Corte 


Anteriormente ya se presentó la metodología de análisis de esta relación, 
la función de CEINCI-LAB que permite este análisis es 
conexion_viga_columna_V2 y dibujo_Capacidad_Viga_Columna representa 
gráficamente la relación en los elementos (ver figura 52a). Como se puede 
apreciar en la figura 52b, todas las relaciones son menores a 1 excepción del 
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último piso, sin embargo, las normativas indican que en este último nivel no se 
requiere cumplir ese requisito. 


Código del incipal, parte 13: 


%% Verificación de la Conexión Viga - Columna 
[sum_M_col,sum_M_vg,Capacidad_Vg_Col]=conexion_viga_columna_V2(Fy,Fu 


,coseno,nr,nudcol,X,Y,NI,NJ,L,Z_ ELEM,Seccion,Elem_Tipo_Prop2,Sec_VG_COL, 
Fuerzas_Elementos CM,Fuerzas Elementos CV, Envolvente_omega); 


dibujo Capacidad Viga _Columna(X,Y.NI,NJ,Capacidad_Vg_Col) 


(a) 





Verificación de Capacidad *Viga/Columna 


Vg/Col = 0.64 | + , g/Col = 0.64 


g/Col = 0.43 g/Col = 0.76 g/Col = 0.76 g/Col = 0.43 


g/Col = 0.47 g/Col = 0.83 g/Col = 0.83 g/Col = 0.47 





g/Col = 0.52 g/Col = 0.94 g/Col = 0.94 g/Col = 0.52 


0.0 0.9 0.95 


(b) 
Figura 53 a) Código de CEINCI-LAB; b) Relaciones Viga. 


6. COMENTARIOS Y CONCLUSIONES 


La filosofía del sistema de computación CEINCI-LAB es desarrollar 
programas en Matlab con el propósito de que los usuarios aprendan como se 
realiza un determinado proceso y concomitante con esto presentar el marco 
teórico que facilite la autoeducación. 


Con este antecedente, en este artículo se han desarrollado nuevas 
funciones de análisis y de diseño de una estructura de acero conformada por 
vigas y columnas, con conexiones con capacidad a momento. Se inicia el estudio 
desde el prediseño de los elementos; posteriormente se realiza el análisis sísmico 
de dos formas: la primera empleando el Método Espectral y la segunda con el 
Método Estático. 
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Una vez que se han determinado las fuerzas estáticas equivalentes 
debido a la acción sísmica, se pasa al análisis estático, considerando 9 
combinaciones de carga, se resuelve cada una de los estados de carga y se 
determina una envolvente donde se seleccionan los valores máximos y mínimos 
de fuerza axial, corte y momento en cada uno de los elementos. El uso de las 
nuevas funciones de CEINCI-LAB , en este artículo a más de presentar el marco 
teórico en forma detallada se indica el uso de las nuevas funciones, mediante el 
análisis de una estructura de acero de 4 pisos y 3 vanos. 


Una vez que se ha determinado la demanda en los elementos, se pasa a 
obtener la capacidad, para ello se presenta en forma detallada la teoría que 
conduce al análisis a flexión, corte y carga axial de vigas IPE y columnas HEB. 
Después se obtienen relaciones demanda/capacidad de cada uno de los 
elementos y se realiza el control de columna fuerte viga débil. 


El uso de las nuevas funciones de CEINCI-LAB hacen que la entrada de 
datos y la visualización de los resultados sea muy fácil de usar y permiten ver si 
un diseño está bien realizado. De tal manera que el lector va a contar con nuevas 
herramientas que le van a facilitar el cálculo estructural. Se aspira de esta forma 
haber aportado al diseño de estructuras de acero. 
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